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Abstract 


The impulsive noise associated with helicopter flight due to Blade- Vortex Inter- 
action, sometimes called blade slap is analyzed especially for the case of a close 
encounter of the blade-tip vortex with a following blade. Three parts of the phe- 
nomena are considered: the tip-vortex structure generated by the rotating blade, 
the unsteady pressure produced on the following blade during the interaction, and 
the acoustic radiation due to the unsteady pressure field. 

To simplify the problem, we confine our analysis to the situation where the 
vortex is aligned parallel to the blade span in which case the maximum acoustic 
pressure results. The 2-dimensional incompressible flow is assumed with uniform 
motion of the blade. The tip-vortex is modelled so that the circulation near the 
tip is rolled into a concentrated vortex and the extreme case of the interaction 
is studied when the following blade cuts through the center of this vortex core, 
which is turbulent and viscous. It is further assumed that, during the interaction, 
there is no distortion of the vortex path or of the vortex itself ; in other words 
the interaction occurs only through the boundary condition on the blade giving 
an unsteady pressure on the blade surface, the usual assumption for gust analysis. 

Acoustic radiation due to the interaction is analyzed in space-fixed coordinates 
and in the time domain with the unsteady pressure on the blade surface as the 
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source of chordwise compact, but spanwise non-compact radiation. Maximum 
acoustic pressure is related to the vortex core size and Reynolds number which 
are in turn functions of the blade-tip aerodynamic parameters. 

Finally noise reduction and performance are considered. Some reduction in 
noise may be achieved by modification of the blade-tip shape but at the expense 
of reduced aerodynamic efficiency. 
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Chapter 1 
Introduction 


The tip vortex for many years has been the subject of research and is important 
to the generation of the lift and drag of the wing and the downwash at the tail 
of an aircraft. In addition research on the tip vortex has been stimulated by 
the concern that the tip vortex of a large aircraft might be hazardous to smaller 
following aircraft. 

Recently, interest has been increased in the tip vortex produced by a rotating 
blade not onl because its path and core size are related to the prediction of the 
rotor blade performance but also because the vortex encounters the following blade 
or the tail blade periodically, and causes the one of the primary sources of vibration 
and impulsive noise, sometimes called “ blade slap”, in certain flight conditions 
for helicopters. 

During the interaction with the blade, there will be fluctuation of the flow 
around the blade and generation of acoustic waves radiating to the far field. It is 
also expected that the strength and the structure of the vortex become important 
factors during the interaction especially for the close interaction with the blade. 

The objective of this thesis is the calculation of the unsteady loads and the 
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acoustic field, arising from the interaction, with the rotor vortex described in 
terms of tip balde aerodynamic parameters, which are directly related to the per- 
formacne, blade tip design and noise reduction. 

To analyse this blade vortex interaction, three parts of the phenomena must be 
considered; the tip vortex generated by the blade, the unsteady pressure produced 
on the following blade during the interaction, and the acoustic radiation due to 
this unsteady pressure field. The vortex is concentrated in a single curved filament 
whose core is viscous and turbulent, and during the interaction with the following 
blade, the vortex and vortex path are distorted and sometimes a secondary vortex 
is generated at the surface of the blade. Meanwhile the pressure on the surface 
becomes unsteady due to the change of local instantaneous angle of attack, the 
generation of the trailing vortices and the acceleration of the flow. This disturbed 
flow and fluctuating pressure on the surface generate the source of the acoustic 
wave propagating to the far field. 

Three extreme case of blade vortex interaction can be considered as shown in 
Figure 1 ; first, the vortex filament is parallel to the blade span, secondly ; the 
vortex is perpendicular to the span such that the filament direction is parallel to 
the blade chord ; and finally, the filament direction is perpendicular to both the 
blade chord and span, as happens during a main rotor vortex and tail rotor blade 
interaction. Only the first case of parallel interaction is analyzed in this thesis. 

In the next section, a general background is presented briefly for the flow 
interaction, aeroacoustics, tip vortex and helicopter noise. Previous work on the 
blade vortex interaction are reviewed in the following section. The problem and 
the approach of this thesis are described in the final section in this Chapter. 
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In Chapter 2, from the general aeroacoustic equation for a moving body, linear 
parts of the equation ( i.e. monopole and dipole sources) are reformulated including 
the effect of a moving observer using the kinematic relation between time and space 
variables. 

In Chapter 3, a detailed analysis is given for the rotor tip vortex, the unsteady 
loads during the interaction and the acoustic pressure due to the unsteady load. 
The noise-performance trade off is also discussed. 

In the subsequent chapter, some results are presented and discussed for the 
vortex velocity profile, the unsteady lift and acoustic signal, and the effect of the 
• vortex structure and observer position. Noise reduction for various tip circulation 
shapes is included. 

Conclusions arising from the results are summarized in Chapter 5. 

Finally, several suggestions for further research on the blade vortex interac- 
tion are made particularly regarding the unsteady compressiblity, noncompact- 
ness, nonlinearity and vortex generation. 
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1.1 


Background 

1.1.1 Flow Body Interaction 

The vortex moving toward the blade can be considered as a convecting disturbance 
to a body. Generally, the convecting disturbances in the flow cause fluctuations 
of the flow around the body and at the body surface during the interaction, and 
this unsteady fluctuating flow causes unsteady forces on the body surface and the 
generation of acoustic waves propagating to the far field. 

This phenomena occurs in a variety of situations in flight vehicles : ; aircraft in 
an atmosheric gust, manuvering missile configurations and canard wing aircraft, 
fans and compressors in a turbo engine or counter rotating propellers, as well as 
in helicopters flight. 

The' convecting disturbance may be turbulence in the atmosphere or wind 
tunnel, turbulence in the jet or wake, oscillatory flow in the mixing layer, a single 
concentrated vortex such as a wing tip vortex or starting vortex, or an incident 
acoustics wave or blast wave. The main flow may be uniform or nonuniform. 
And the body may be an airfoil, plate, edge, corner or plane normal to the flow 
direction. 

As summarized by Rockwell^!, during the interaction, several kinds of complex 
phenomena can be expected , which are difficult to analyze : secondary vortices or 
instability waves may be generated, laminar - turbulent transition may be triggered 
and the incident acoustic wave may amplify these fluctuations. 

However, for the case of a small disturbance due to the interaction, the velocity 
field may be split into the vorticity field representing the oncoming disturbance 
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and the potential field representing the disturbed field^’^l. The potential field is 
directly related to the pressure field on the body surface and to the acoustic wave 
at the far field. The vorticity field may be distorted or not during the interaction 
depending on the flow and body conditions. Whether distorted or not, the two 
fields are coupled at the body surface to satisfy the impermeable conditon of the 
body surface, which condition makes the streamlines of the flow distorted^! as 
shown in FigureZ. The above mentioned concept has been used in conventional 
gust theory!^ and is analogous to the rapid distortion theory of turbulence^. 

1.1.2 Aeroacoustics 

For the acoustic field associated with the flow-body interaction, the above the- 
ory can be used^l and extended for strong interactions. Alternatively Lighthill’s 
acoustic analogy approach^ ^ can be used ; this was originally developed to de- 
scribe the acoustics generated by a fluctuating fluid flow occupying a limited part 
of a very large volume of fluid of which the remainder is at rest. The problem of 
jet noise^l represented by a quadrupole source in uniform acoustic medium at 
rest was treated in this way by Light hill. 

Curled extended the analogy for the presence of the body fixed in the flow or 
vibrating about a mean surface by using Kirchhoff ’s formula. He showed that the 
body effect, including the reflexion and diffraction at the surface, can be incorpo- 
rated into dipole sources and these dipole sources become increasingly important 
for flows of low Mach Number. 

Ffowcs- Williams and Hawkings^l generalized the theory for the moving body 
in an acoustic medium at rest represented by moving monopole sources in addition 
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to the dipole and quadrupole sources. For the monopole and dipole sources, his 
theory is the generalization of Lowson’s acoustic theory! for moving acoustic 
singularities. 

For the problem of body fixed in the uniform flow, a Lorentz-type of 
transformation! 1 **! can be used to obtain the acoustic field as shown by 
Blokhintsevl^l. For the varying acoustic medium, as for a shear flow, the acoustic 
analogy was reformulated by Phillips! 1 **! and Lilley! 1 ®! to account for the source 
convection and the sound refraction. 

1.1.3 Tip Vortex 

The tip vortex, as the result of previous research, may be described as follows : 
Streamwise vorticity sheds from the trailing edge of the wing to form a continuous 
vortex sheet which rolls up into a concentrated tip vortex and this tip vortex 
remains for long distances behind the wing until it decays and finally becomes 
unstable as shown in Figure 2. The concentrated vortex has a viscous region near 
the center of the core, a turbulent region for most of the core and potential region 
outside of the core. 

The assumption that the vortex sheet remains flat has been generally accepted 
in the calculation of the lift of high aspect ratio wing, t.e, the well known lifting 
line theory of Prandtl. However, in the calculation of the downwash at the tail 
of an aircraft!^!, the assumption is too restricted especially for low aspect ratio 
aircraft because the sheet rolls up quickly within a chord length from the trailing 
edge. So researchers has been interested in the rate of roll up [22], [25] an( j the 
core dimensions of the vortex. It is known that the vortex sheet rolls up quickly 
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for a low aspect ratio wing at large angle of attack^, and that the circulation 
of the vortex increases to the maximum circulation at the root of the wing if the 
vorticity shedding from the wing is assumed completely roll up. 

Prandtll 2 ®! calculated the vortex core size with the assumption of solid rotation 
in the core and potential flow outside of the core, using that the relation between 
rotational kinetic energy in the wake and the induced drag. Based on the calcula- 
tion, the core size is 0.086 of a span length] 2 ®! for the elliptic loading distribution. 
Spreiter and SacksI 2 ^ used the same assumption and caculated more carefully 
and showed that the size is around 0.078 of the span length. Due to the simple 
model of the uniform vorticity core, Donaldson! 2 ^ pointed out that the core size 
is overestimated and the peak velocity is underestimated. And he showed that 
the Betz theory! 21 ! predicts the overall velocity profile of the vortex reasonably 
well. Betz developed the theory with the assumptions of conservation of vorticity, 
angular momentum and second moment of inertia in the circulation between the 
wing and the vortex. However due to the inviscid theory of Betz, the velocity at 
the center becomes singular for elliptic loading and there is an ambiguity in the 
core size in general. 

With the advent of large aircraft, the research! 2 **! on a more detailed vortex core 
structure and the' decay of the vortex became important and Squire! 2 ^! obtained 
a solution for the turbulent viscous vortex with the assumption of constant eddy 
viscosity by using the Lamb’3 solution^! for the laminar core. He showed that 
the turbulent core decays faster than the laminar vortex. For the structure of the 
vortex, Hoffman and Joubert^^ obtained the log law profile of the circulation 
in the inner region of the vortex by applying the mixing length theory to the 
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vortex. The profile was also confirmed with experiment data. Saffman^l’l^l, an( j 
Govindarju and Saffman^! tried to explain the general properties of the turbulent 
vortex core without a specific turbulence model ; they characterized the vortex 
in terms of possible overshoot in circulation profile, rate of decay of maximum 
velocity and slow decreasing function of Reynolds Number. They^ also pointed 
out that some data measured for the maximum circulation of the vortex is only 
half of the circulation at the root of the wing. Detailed properties of the turbulent 
vortex were shown by Iversen^L He calculated the angular momentum equation 
numerically with a variable eddy viscosity along the radial distance and showed 
the constant behavior of the core and no overshoot at large Reynolds Number. 
Philips^®! estimated the Reynolds stress based on the velocity profile obtained by 
matching the log law of Hoffman and Joubert near the peak velocity and inviscid 
profile at outside of the core obtained by Moore and Saffman^l for the initial 
stage of roll up. 

In mostof the above analysis, the axial velocity is not considered, however, it is 
believed that it becomes important in instability of the vortex and in calculation 
of the drag. Bachelor^ calculated the drag due to the axial velocity defect as 
well as due to the rotational vortex for the laminar case with a solid rotation core 
modeh. Recently, Roberts^®] showed that the axial velocity has an influence on 
the initial persistence of the concentrated vortex and on its subsequent decay. He 
explained that this is due to the balance of the turbulent diffusion by a radial 
inflow caused by the acceleration of the flow in the core. He also obtained the 
velocity profile by solving the turbulent diffusion equation in closed form. 
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1.1.4 Helicopter Noise 

Helicopter noise, together with propeller noise, has been recognized as important 
because of the environment problem during flight close to residential areas or in 
cities and because of the interior noise in civil helicopters. In addition, helicopter 
noise is significant in the eariler detection during day and night operations of 
military helicopters. So the understanding of the noise source mechanism and the 
determination of parameters which can reduce noise and influence on the blade 
design within the restriction of good performance are of current urgent interest. 

There are two kinds of helicopter noise in general ; discrete noise and broad 
band noise in the frequency spectrum as summarized by LowsonI 4 ®! and later 
by George^ 4 !. The discrete noise, having fundamental frequencies according to 
the blade passage frequency and its harmonics as shown in Figure 4, is due to 
the rotational steady loadsl 44 !, azimuthly varying loadsl 4 ^, volume displacement 
of airl 4 ^l by the high speed rotaing bladet 4 ^’^ and blade vortex interactions. 
Especially, the noise due to the thickness, and the blade vortex interaction, called 
the impulsive noised, have a large number of harmonics and the signals are 
observed impulsively in time as shown in Figure 5. The thickness noise is mainly 
due to the rotation of the blade causing the propagation of a crescent shaped 
wave or local shock from the rotating blade tip^l to the far field for a high speed 
rotor. The blade-vortex interaction noise is mainly due to the rapid load variation 
during the interaction. The broad band noise, observed between the harmonics 
in discrete noise and at the high frequency band above around 500 Hz, is due 
to the turbulence in the boundary layers, the vortex shedding at the tip and the 
inflow of the atmospheric turbulence causing the high frequency trailing noise^l 
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and low frequecy interaction noised. Both the discrete and broad band noise 
are important. However, considering the threshold of hearing and atmospheric 
attenuation for far distances^, the band width in the range 250 to 500 Hz is the 
most annoying. 
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1.2 Review of Previous Work 


This blade vortex interaction phenomena was first studied by Leverton and 
Taylor^l experimentally using a rotating blade and two opposed airjets to simu- 
late the tangential velocity profile of the vortex. Using the theory of nonuniform 
flow past a thin airfoil formulated by Von Karman and Sears^i, Sears^l treated 
the line vortex as a gust where the vortex is forced to move parallel to the blade 
span with a displacement of half of the chord height ; Sears calculated the acoustic 
signal by using a dipole source to represent the unsteady force on the surface. 

The restriction to the forced vortex was relaxed by Parthasarathy and 
Karamcheti^J to analyze the unsteady field due to the interaction between the 
blade and a free vortex, whose path is disturbed by the blade during the inter- 
action from the initial height of half chord, and they calculated the quadrupole 
effect corresponding to the acoustic sources in the flow field. 

Widnall formulated an acoustic model for an oblique forced vortex using the 
quasi 2-D unsteady aerodynamics for an oblique sinusoidal gust in incompressible 
and subsonic compressible flow obtained by Johnson^l and Filatos 1^1. The 
Fourier transform of the velocity of the vortex and the pressure field was used 
to analyze the acoustic farfield and the viscous core in the potential vortex was 
treated by using an effective distance between the center of the vortex and the 
blade. Further, the Betz inviscid vortex model^^ was used^J to relate the vortex 
core to the tip loading shape. Filatos^ also studies the blade loads and the noise 
due to interactions with an array of equally spaced oblique vortices. 

For the high speed rotor, noncompactness becomes important especially for 
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high frequency tones, where the acoustic wave length becomes smaller than the 
characteristic body length i.e. chord length. For this reason, high frequency ap- 
proximate unsteady subsonic theories for the sinusoidal gust were developed by 
Adamczyk! 7 *! and Amiet!®®! separating the problems at the leading edge and the 
trailing edge as proposed by Landahl!®®!. The theories were accurate for gust 
wave length A < times chord length, when compared to the exact solution 

obtained by Graham!®®! using the similarity rules ; these were extended by Mar- 
tinez and Widnalll^l’^l, including the acoustic field, for subsonic trace speed 
M tr = Moo / sin A < 1, where A is the oblique angle of the vortex. Kerschen!®^! ex- 
tended the acoustic theory for A/< f > 1 and studied the steady loading effect!®**! to 
the unsteady fields by using the rapid distortion theory generalized by Hunt!®! and 
Goldstein^ 7 !. To study the vortex tail rotor interaction, Amiet!®®! developed the 
unsteady response function for a delta function gust, which was done by Meyer!®^! 
for incompressible flow, and explained the effect of noncompactness on the acoustic 
field by using the dipole source. 

With the recent development of the numerical schemes, a number of researchers 
have studied the nonlinear effect due to the the vortex induced shock at transonic 
speeds!®®!’!®®!’!®®! anc j due distortion of the vorticity field at low Reynodls 

Number!® *1. Others!®^! have tried to predict the noise of a full scale rotating blade 
using the Ffowcs- Williams and Hawkings equation!^! with the measured unsteady 
aerodynamic pressure on the surface!®®!. Several full scale measurements of the 
noise, measured at the helicopter!®®! or other aircraft moving with the same speed 
of the helicopter!®'*! or a fixed location of the observer!®®!’!®®!, have been reported. 
Small scale rotating blade experiments!® 7 ! from the acoustic wind tunnel were also 
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conducted and the signals were compared with those in full scale measurement. 
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1.3 Statement of the Problem 


Prediction of the unsteady pressure and noise generated by the blade vortex in- 
teraction in full scale is beyond from accurate theoretical method^ 2 !. Even ex- 
periments using small scale models are difficult due to the lack of correct scaling 
parameters. This is because the phenomena depend on the flight operating con- 
ditions and blade geometry which in turn depend on the vortex strength and 
structure, its position before the interaction, and its path during the interaction 
as well as the Mach Number and Reynolds Number of the rotating blade. In con- 
trast to the fixed wing, research on the rotor vortex structure^ is less complete 
and quantitative datal 4 ^’! 42 !’! 4 ^ is very limited. 

In previous anaysis, the above mentioned works used a point vortex or an 
inviscid vortex which has a singularity at the center of the vortex core. Since the 
vortex strength and structure become important for the case of a close encounter 
with the following blade^l j n certain flight condition as shown in Figure 7, a 
more realistic vortex should be used. 

It is attempted here to study the effect of the turbulent viscous vortex core on 
the unsteady loading and the acoustic pressure field. The vortex is described in 
terms of the Reynolds Number, T/v and the turbulent vortex size , r 0 , where T is 
the constant circulation out side of. the vortex core. To complete this phenomena, 
the vortex strength and structure are determined from the blade tip aerodynamic 
parameters. 

To simplify the problem, the situation is considered where the vortex filament 
generated from the rotating tip is aligned parallel to the following advancing blade 
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as shown in Figure 7 ; in that situation the maximum acoustic pressure due to 
the interaction is observed 1^3] because of the large interaction area and in-phase 
radiation^!. Even with the assumption of 2-D aerodynamics, the flow is too 
complicated to analyze completely as indicated in the experiment done by Ziada 
and Rockwell^! and recently by Yu^^aand Mandella, Moon and Bershaderl^l 
; thus it is assumed that during the interaction the vortex path and the vorticity 
field are not distorted when the vortex is cut by the following blade. The blade 
is assumed thin compared to the turbulent vortex size and thick enough to avoid 
separation at the leading edge. This is the linear assumption used in conventional 
gust theory and rapid distortion theory of turbulence. The flow is assumed in- 
compressible for the interaction in the middle of the span as shown in Figure 7. 
Even though compressible effects become important for the high frequency region, 
there is no significant effect of the compressibility as shown by Graham^, for the 
disturbance wave length A greater than one and half of the chord length at Mach 
Number around 0.4. 

Acoustic pressure is calculated in space fixed coordinates and reception time 
domain, whereas the source is described in the body fixed coordinate and emission 
time domain. The source is assumed chordwisely compact for the above conditions 
where the significant disturbance wave length is larger than the chord length at 
moderate Mach Number around 0.4. (Consider that the observed acoustic wave 
frequency containing most of the energy due to the blade vortex interaction is from 
zero to 750 Hz^’^1 where the most significant acoustic wave length is around 
lm and that the most annoying frequency is around 250 to 500 Hz as shown in 
Figure 6.) However, the source is assumed spanwisely finite i.e., noncompact, 
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where the retarded time differences exist between elements along the span length'; 
Thus the noncompact source will influence effects the signal for an observer either 
fixed or moving with the same velocity of the source. 

Based on the analysis discussed above, the trade-off between noise and perfor- 
mance is determined by varying the shape of the blade tip circulation. 
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Chapter 2 

Aeroacoustic Formulation 


The acoustic analogy approach to determining the sound generation due to the 
fluid motion in an atmosphere at rest was introduced by Lighthill. The induced 
acoustic field is determined from the prescribed nonlinear flow as a quadrupole 
source in the wave equation. This approach was extended by Curie to include 
the body surface effect. For a surface in arbitrary motion, Ffowcs Williams and 
Hawkings^ ^1 derived the equation by using generalized functions under the as- 
sumption of an impermeable surface. The solution shows the Doppler factor due 
to the surface motion derived in the time domain. The Doppler shift in frequency 
due to the source motion is not shown explicitly. However the time and space 
relation used in the solution implicitly contains the Doppler shift. Moreover, the 
Doppler shift due to the observer motion, which was not discussed in the above 
approaches, can be derived from the time and space relation with the source and 
observer in arbitrary motion. In this chapter, the Ffowcs- Williams and Hawkings 
aeroacoustic equation will be discussed with the more general surface conditions. 
The effect of both the source and observer motion will be discussed more specif- 
ically with alternative expressions for the thickness term and the dipole term. 
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These are additional to the quadrupole term and are due to the surface motion . 


2.1 General procedure for the aeroacoustic 
equation and solution 


The general compressible and viscous flow motion can be expressed in the Navier- 
Stokes equation as follows. 

Continuity equation 


dp dpUj 
dt dxi 


m( x it t ) 


( 2 . 1 ) 


Momentum equation 


dpui , dpmuj di-pSij + n,) , ^ 

-W + SZ, + ' ,/ ' (z " i) 


( 2 . 2 ) 


where 

p density of the fluid 

u, velocity of the fluid in i direction 

p pressure of the fluid 

m mass change rate in the fluid element 

/ body force in the fluid element 

8{j Kronecker delta (S,y = 0, i ^ j = 1, * = j) 

Tij viscous stress tensor tv/ = p. (gi + - ffygj) 
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The surface in motion should satisfy the following equation 


d ± 

dt 



where the surface is defined by the equation 


(2.3) 


/(*(,<)= 0 (2.4) 

and V{ is the velocity of the surface. 

The above continuity and momentum equations are valid in the fluid, that is, 
outside of the surface /(z,-,t) = 0. Thus, if there is a body, the equations outside 
of the body surface can be expressed by using the Heaviside function defined as 


H[f) = 0 

/< o 

(2.5) 

= 1 

/> 0 

(2.6) 


where /(z,-, t) > 0 is for the outside of the surface and /(z,-,t) < 0 is for the inside 
of the surface. 

The interesting quantity here is the perturbed density ( p — p 0 )H(f) where p 0 
is the density in the medium. The continuity equation (2.1) becomes as written 
below for the region outside of the surface, after multiplying equation (2.1) by 

*(/), 


§1 {{P~Po)H } + 


= Po6tf)vi^ + p6(f){u { -Vi)^ + m(xi,t)H (2.7) 
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where 5(f) is the delta function having zero value except at / = 0, that is, at the 
surface. 

With the same procedure, the momentum equation (2.2) becomes 
d d 

— + — {(pttftty + p6ij - Tij)H) 

= {/»«.■(«; - v j) + Phj ~ + Pf' H ( 2 - 8 ) 

The above equations can be obtained by using the properties of the Heaviside 
function described in Appendix A.l. and by using the equation of the surface 
motion (2.3). 

From the two equations (2.7) and (2.8), the aeroacoustic equation can be ob- 
tained by the following procedures : 

Take the operator ^ for the continuity equation (2.7) and the divergence 
for the momentum equation (2.8). Add the two equations and substract in 
both sides of the combined equation. Then the aeroacoustic wave equation is 



where 

Tij = pu^ + {p - a 2 0 (p - p 0 ))5ij - Tij (2.10) 
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Fa 
Qi 

Tij is the Lighthill quadrupole strength and F,y is the dipole strength and is 
the monopole strength. 

The above inhomogeneous wave equation is nonlinear and there is a term in- 
cluding the dependent variable {p - p 0 ) in the right hand side. However if the 
right hand side is considered as the source term prescribed beforehand , the wave 
equation behaves like a linear inhomogeneous equation with the sources in motion. 
The solution of the equation (2.9) can then be easily obtained as following : 

= lJ r jg^.mfWAyi,T))H9)dy<dT 

+ f, l TiJr T ) + %)<%*• 

(2.13) 

where R is the relative distance between the source and the observer 

. R =| i, - y, | (2.14) 

and g is the function describing the space and time relation as below 

R 

g = t — t H (2.15) 

&0 


= -{phi ~ Tij + ~ V;)} (2.11) 

= poVi + p{ui - Vi) ( 2 . 12 ) 
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The coordinate Xi,t is the position and time of the observer and the coordinate 
y,-, r is the position and time of the source in the space fixed coordinate as shown 
in Figure 8, where the observer and source may be in motion or fixed. The 
coordinates have the following relationship : 


d 

dyi 

d 

dx{ 

(2.16) 

d 

d 

(2.17) 

dr 

dt 


By using the above relation and by taking the partial integration over space y, 
for quadrupole and dipole terms and time r, the solution (2.13) becomes 


4*- a 2 0 H(f)(p- p 0 )(xi,t) = 


+ 

+ 


d 2 

dxidxj 


j l^mf)Td*.r)M9)*yidT 


17,1,1,1 (W) r > - o) % )iydr 

Itl, ll T * + r )) Hn)dy.dr 

(2.18) 


In the solution , the sources are defined in the space fixed coordinate y,-. How- 
ever it is useful to specify the sources in the body fixed coordinate r/,- moving with 
the velocity v as shown in Figure 8 such as 


y.- = y, + J Vi(n, t'W 


(2.19) 


The volume changes of fluid element in the two coordinate system gives 
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dy { = diji 


( 2 . 20 ) 


The volume integration for the dipole term including Fij and the monopole 
term including Q, in the above equation can be changed into a surface integration 
and the time integration over r in the solution can be carried out by using the 
properties of the generalized function described in Appendix A. 2. 

From equation (2.15) , (2.14) and (2.19), g has the following form for both the 
source and observer in arbitrary motion 


R 

g = r - H (2.21) 

d 0 

where 

R = | - y,- 1 

= | x ° + J v i{xi,t')dt' — rji — J v t (r')dr' | (2.22) 

v is the velocity of the body and v° is the velocity of the observer and x° is the 
initial position of the observer. 

Then 


dg , dR 1 

= 1 + 

or ot a 0 

t Xj - y,- t;,- 

I Xi - y,- | a 0 

= 1 -M r 


(2.23) 


where M r is the source Mach component in the direction to the observer as shown 
in Figure 8. The above relation depends on the source velocity v,- explicitly 


23 


whereas the observer velocity v° is implicitly related in the relative distance | 
X{ — yi | between the source and the observer. 

Thus, the final form of the solution (2.13) becomes as given below by using the 
equations (A.6),(A.8),(2.20) and (2.23) 


*(/)(/>-/>.)(*»*) = 


a 2 ! ftifo.r.) 

dxidx.- Jv F 4ira 0 R I 1 — M r 

4 i r c 


_9_ f [ PfijVii T e ) 

dxi Jv F [47ra 0 i2 | 1 — M r |_ ^ ' 

d_ f [ m(r?,-,r e ) 

dt Jv F [4ira 0 R | 1 — M r | u * 

_d_ f [ ds 

dii Js [4ira 0 ^ | 1 - M r | Tc 

f UQi[Vi,Te) 

dt Js 4ira 0 R \ 1 — M r | r 


where /,• is the outward normal direction to the surface S defined as 


(2.24) 

(2.25) 

(2.26) 

(2.27) 

(2.28) 


df/drjj 
I df/drn | 


(2.29) 


V F is the volume actually occupied by the fluid, 5 is the surface of the body and 
the integrands are evaluated in the body fixed coordinate rji and at the emission 
time r t defined in (A.9) and (2.21). 

The above expression for the acoustic solution is the general form including the 
quadrupole term, T, ; - and the dipole term, and the monopole term, sometimes 
called thickness term, and terms of mass change and body force. The surface is 
not necessarily impermeable, that is, which normal velocity /,v,- is not necessarily 
the same as the normal velocity of the fluid The body and the observer can 
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be in arbitrary motion. However, this is restricted to a homogeneous medium at 
rest. 

For the case of the moving surface considered here, the monopole and dipole 
terms dominate the quadrupole term at the relatively low Mach Number. The 
source strength Fy and Qi in (2.11) and (2.12) can be simplified for the rigid surface 
such that (u,- — v t ) = 0. In the following section, the alternative useful expressions 
for the monopole and dipole terms for the rigid surface will be discussed in the 
case of no mass change and no body force in the fluid. The kinematic relationship 
between the time and space for both the source and observer in arbitrary motion 
will be also discussed. 

2.2 Kinematic Relationship between the Space 
and Time 

From equation (A.9) and (2.21), the time and space relation between t, r e and R 
becomes 

T t — l 

a 0 

where 

R =\ x° + J v D (xi, t')dt' — rji — J r')dr' | (2.31) 

Then the relation between the time scales in emission time r e and in reception 
time t can be obtained from the following 
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dR 1 


Thus 


dr e 

dt 


dt a 0 

R±Vi_ Ri Vi dr e 

R a 0 R a 0 dt 


dr t _ 1 - Rj/R ■ v°/a 0 

dt \-RijR-Vija 0 
1 -M° r 
1 -M r 


(2.32) 


(2.33) 


The relationship between the time scales of the equation (2.33) can be related 
to the conventional Doppler shift in frequency as follows for the special cases that 
the motions of the source and the observer is in the same straight line or the 
velocities of the source and the observer is the same. 

In that cases, ^ is constant, so 


r e 


1 - M° 
1 -M r 


t + Ci 


(2.34) 


where C j is constant. 

For a simple harmonic source, the source is described in the emission time 
domain r e as • 


h = he iu ' r ' (2.35) 

where u e is the source frequency. 

This source will be observed in the reception time domain t in the form below 
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p = \pe VJtTe \ Tc =t-R /ag 


(2.36) 


By using the equation (2.34), the pressure p will be 


- * (“ 
p = pe v 


1 zOjL t ur r q I 

l-A/ r 1 — A/ oo / 


(2.37) 


From the above equation, the observed frequency u a becomes 




1 - M r ° 
1 -M r 




(2.38) 


which is the conventional Doppler shift in frequency. 

The above relation in frequencies comes from the kinematic relation in time 
and space (2.30). ( The relation can be also obtained from Figure 9 for the special 
case of linear motion). It can be said that both the observer and source motion has 
an effect in frequecy , but only the source motion has an effect on the amplitude 
through the factor | 1 — M r | in equation (2.23) and (2.24). In other words, the 
observer in motion just collects the wave information already accumulated due to 
the source in motion but differently from the stationary observer in frequency. 

2.3 Alternative expressions for the monopole 
and the dipole terms 

2.3.1 Monopole term 

The monopole term in equation (2.13) becomes as shown below by using (A. 3) 
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= f 

dt J yi Jr 

= ~t I 

dt J Vi Jr 


po—^~b{g)dyidT 


4ira 2 R dr 


^^{1 - H(f)}S( g )d yi dr (2.39) 


Hereafter, ( p — p 0 ), itself, represents the density fluctuation in the outside of 
the surface. 

Using equations (A. 8) and (2.23) and taking the partial integration over the 
time r as done in (2.18), equation (2.39) becomes 

( P~Po){Xi,t ) 


where V B is the volume occupied by the body, that is, the volume inside the 
surface. 

The above expression is useful especially for the case of compact source. In 
that case, equation (2.40) becomes 


-( [ 

dt 2 4 Jr 

d*_ C 

dt 2 Jv B 


Mi - #(/)} 

4iralR 

Po{l~H(f)} 


6{g)dyidT 


Ava 2 0 R 


1 -M r 


drji 


(2.40) 


Jr. 


d 2 

( P~Po){Xi,t ) = — 


PoV B 


AiralR | 1 


Mr 


j r. 


. (2-41) 


which shows that the strength of the monopole is the same as the displaced air 
mass p 0 Vb having the same volume of the body. 

From the original expression in (2.24) it is hard to describe the compact source 
strength because the strength of the monopole in (2.24) is expressed in terms of 
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the normal velocity of the surface having the positive value in the front of the body 
and the negative value in the rear of the body for the closed surface. This situation 
occurs also in describing the thickness effect in aerodynamics. The thickness effect 
can be represented by the distribution of sources and sinks at the surface of the 
body ; or it can be represented by a doublet distributional for which the direction 
vector is opposite to the direction of the uniform flow for an axisymmeric body, 
or a single doublet for which the strength is related to the area for a cylindrical 
body and related to the volume for a sphere. 

Unlike the stationary monopole source, which has no directional pattern for 
a compact source, the moving monopole source has a directional pattern in the 
direction of the source. This means that the source representing the effect of the 
thickness in motion has the property of a dipole source whose direction is the same 
as the direction of motion, as in the case of the aerodynamics explained above. 

The equation (2.41) was used in the calculation of the acoustic pressure due 
to the rotaing body in the paper by Wright and Lee^l eV en though the equation 
was derived in different wayl^l. It turned out that there is no difference between 
the chordwise compact source and the noncompact source up to Mach Number^ 
0.7. 

The more useful expression can be obtained by changing the time derivative in 
t into r e having the same variable inside [ ) in equation (2.40) or (2.41) such as 


( P-po){Zi,t ) 


dr e d j 

f dr e d 

1 Po^B 

J 

dt dr e ' 

\ dt dr e 

4tt a 2 0 R | 1 — M r |_ 


(2.42) 


Using the equation (2.33) for 4^-, the density in equation (2.42) can be calcu- 
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lated in the emission time r e domain as 


(/>-/>•)(*<■ o 


1 -M f ° d 

( 1 — M° d 

1 PoVb 

J 

1 - M r dT e 

{ 1 - M r dr e 

4ira^R | 1 — M r 


(2.43) 


The density in the observer time t domain can be transformed by using the 
time and space relation of the equation (2.30). In the case of the noncompact 
source, the density can be also calculated from the equation (2.40) with the same 
procedure used in the case of compact source. 

2.3.2 Dipole term 

The dipole term in the solution (2.13) is 


{p~Po){Xi,t) = 


JL [ l i F n{Vii r <) 
dxi Js 4ira 0 R | 1 — M r 


(2.44) 


For the moving body, it is convenient to change the space derivative in the above 
equation to a time derivative. First, let’s consider the compact source. Then the 
equation becomes 


(P ~ Po){Xi,t) 


d 

dxi [4ir a„R j 1 — M r j 


PM 




(2.45) 


where 


F = IF ■ 


(2.46) 


Using the following relation for arbitrary function F(z,-, T e (x,-, £)) = F(xi,t), 
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F{x u t) I - — F(x { ,T e (x if t)) + +—F(z i ,r e (z i) t )) 


dz{ 


dr t 


(ft), 


(2.47) 


where the subscript t in ( )< indicates that the variable t held constant, equation 
(2.45) becomes 


t P-Po)(x { ,t ) = 


+ 


± E ft(r.) 

dzi 4na 0 R | 1 — M r | 


d_ 


4ita 0 R | 1 — M r | 



(2.48) 


From equation (2.30), 



in the above equation can be obtained as 




dT e A Ri 1 

dzi) t R a 0 \ \ — M r \ 


(2.49) 


The equation (2.45) has two terms : the first term has the space differentiation 
and the second term has the time differentiation. The first term becomes 


_d_ Fj 

dx{ R | 1 — M r 


where 


1 dR Fj Fj d 
R 2 dzi | 1 - M r | R dzi 1 

1 [Ri 1 (Ri RjMj 

R 2 \ R + | 1 - M r | V R R 



(2.50) 

Fj 

1 - Mr | 2 


dR 

dzi 




R<_ 

R 


-±( %^7r L -M i ) 

R V R R J 


(2.51) 

(2.52) 
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The second term becomes 


d_F[x 

dx t R 


where 


i,r(xi,t)) ( dr e \ = 1 d 1 dFj 

1 — A/ r | {dxj, '|1 -M r \ 2 dr c l r + 1 1 — M r | 3r e 

1 , 1 // RM\ 2 >ir2 \l 

= -Fj— + iT^((-R-J -"*JJ 

(2.53) 




1 Ri 


dF { 


R | 1 - M r 2 R R 1 - Af r | 2 dr e 




5R 

5r e 


R, 

-R y *- 


(2.54) 

(2.55) 


Combining (2.50) and (2.53), the equation (2.45) becomes 


{p~ Po)(Xi,t) = 


1 Ri 1 


dFi 


4iral R | 1 — M r | 2 dr e 


+ 0(R~ 2 ) 


(2.56) 


where 0(R~ 2 ) is the second order for the distance R , which can be negligible for 
far field condition and for unsteady force F,(T e ). 

For the noncompact source, it will be 

if - *><*■ « — i h.L miTn ' 7i)l dSM + 0[R ' 2) (2 - 57) 


This final form will be used to calculate the acoustic field due to the blade 
vortex interaction considered here. The unsteady load due to the interaction will 
be analyzed in the following chapter with the tip vortex model determined from 
the blade tip geometry. 
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Chapter 3 
Analysis 

3.1 Rotor Tip Vortex 

The most significant interaction between the rotor tip vortex and the following 
blade occurs in forward and descending flight of the rotorcraft. In that case, it 
is expected that the tip vortex core structure becomes an important factor in 
analyzing the interaction. It is observed that the tip vortex quickly rolls into a 
concentrated vortex and persists for many span lengths before it decays (as in 
the case of the fixed wing tip vortex). However, because of the unsymmetric 
disribution of circulation toward the tip and the unrolled inboard vortex sheets of 
the rotor blade as shown in Figure 10, the rotor vortex is more complicated to 
describe. Thus, for simplicity, the rotor vortex is modelled as an equivalent rotor 
tip vortex pair, i.e., a tip vortex with an equivalent inboard counter vortex of the 
same strength, which replaces the unrolled inboard vortex sheet. The curvature 
of the vortex filament is neglected locally during the formation of the vortex and 
the oncoming velocity near the tip is assumed equal to the velocity at the position 
of the maximum circulation on the span. With this tip vortex model, the vortex 
core structure and strength are related with the tip aerodynamic parameters and 
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are determined by following the general approach used by Spreiter-Sacks^I for 
the potential part and by Roberts!*^ for the turbulent viscous core. Once the tip 
vortex is determined, it is assumed that the vortex persists until it decays. The 
position of the vortex in general depends on the span disribution of circulation . 

3.1.1 Relation between tip aerodynamic parameters and 
vortex core 

Consider a rotor blade of radius R, rotating with an angular velocity fi in a stream 
of uniform velocity U^. The blade will have a point of maximum circulation at 
a radial distance Rm. With the assumptions of an equivalent tip vortex, it is 
convenient to consider the lift and drag on the rotor outboard of the radius Rm. 
Thus, 

L, ip = i,((l R.fCuA, (3.1) 

where A t is the area of tip outboard blade. And the velocity of the blade at the 
radius Rm will depend on the azimuthal angle rp ; i.e., 

U m = (nRm)-(l + ^-smrP) (3.2) 

so that the lift coefficient Cu will in general depend on {U^/URm) sin rp. Alterna- 
tively, the lift can be expressed as the integral of the distribution of tip circulation 
T t around the blade; i.e., 

L,„ = ■ (R - jy r m jf 1 G,d(g~^) (3.3) 
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where T m is the circulation at Rm and G t = T t /T m and the origin of y coordinate 
is located at R ' m as shown in Figure 10. From this alternative expression of lift 
and the equation (3.1), it can be seen that 


I\n = 


Cu 1 
AR t s t 


(3.4) 


where 


• AR ‘= 


(2(R-R m )) 2 2(R-R m f 


2 A t 


At 


(3.5) 


The quantity s t represents the load characteristic of the tip segment of the blade 
and AR t is the equivalent aspect ratio with equivalent span length of 2 (R — Rm)- 
With the assumption of the equivalent tip vortex the drag near the tip can be 
obtained locally from the conservation of the linear momentum in the cylindrical 
control volume encompassing the equivalent blade of the spanlength 2(R — R m ) as 
shown in Figure 12. Neglecting the contributions of the viscous shear stress on 
the control surface, the force on the blade is 


F = —J J pv(v ■ n) dS — j jpndS 


(3.6) 


where v is the flow velocity vector which components are u, v and w in x, y and z 
directions, respectively and p is pressure and n is unit normal vector on the control 
surface S. The x direction force, corresponding to twice the blade tip drag, is 
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D = / J s (Pi~P 2 + P u l 2 - pu 2 )dS 2 - J J s pu 3 v r dS 3 (3.7) 

where subscript ,=1,2,3 represent the quantities on the control surfaces and subscript 
r represents flow quantities on the radial direction of surface S 3 . Si and S 2 are the 
surfaces normal to the x direction, which have same areas, and S 3 is the cylindrical 
surface. In the surface S x , the flow is assumed uniform with Ui = U m , = iwi = 0 
and the surface S 2 is assumed Trefftz plane where w 2 is the twice w at the blade 

i 

trailing edge. For a cylindrical surface of infinite radius, the term involving v r is 
negligible. By introducing the total head H in incompressible flow, the equation 
(3.7), then, becomes as below. 

D = J J H, -Hi+ £(u,’ - u 2 ’) + |(«, 2 + wS)dS 2 (3.8) 

where 

H = p + \p(u 2 + v 2 + w 2 ) (3.9) 

Assuming that the terms involving the differences in H and u between Si and S 2 
planes are mainly due to the viscosity in axial direction x, the last term, which is 
related to the rate of formation of rotational energy in S 2 plane, can be defined as 
the lift-induced tip drag Du. Thus, 

Du = -p(flRm) 2 CDuA t ~ J J s y w 2 2 )dS 2 (3.10) 


where 
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(3.11) 


o<i = p(R - JUr„ / 0 ‘ a,wd( R * R j 

and the induced drag coefficient Con is, in general, a function of (Uoo/QRm) sin 0. 

If the kinetic energy is conserved until the vortex rolls up completely as assumed 
by PrandtlS 0 ! or Spreiter and Sacks^l in the fixed wing case, equation (3.10) 
becomes 


Di = [ f ]-p(v 2 + w 2 )dS + 2 [ \pvg 2 2nrdr (3.12) 

J J pot. 2 Jo 2 


where 


v g = 


2 irr 


(3.13) 


The rotational energy has two parts ; one from the potential part outside of the 
equivalent tip vortices and the other from the vortex cores of radius r„. The vortex 
core structure is assumed not affected by the vortex pair, so it becomes a circular 
flow. The vorticity in the core is assumed a function of the radial distance r rather 
than a constant of solid rotation as assumed by Prandtl or Spreiter and Sacks, 
where the vorticity f = l<9r /dr. The strength of the point vortex, r o , is assumed 
the same as the maximum circulation on the blade, r m , which is actually almost 
90% of r m as determined from the experimental data of Tung^l. 

The surface integral of the potential part can be transformed into a line 
integrals'll by means of the Gauss thorem with the potential function <p, then 



where n is the unit normal vector to the contour of integration following the 
circular core vortex and a branch cut placed between the vortex pair as illustrated 
in Figure 12. It is assumed that the line integral at the far distance is neglected. 

The potential function <j) is considered as a two dimensional vortex pair in the 
Trefftz plane as 


JL tL -1 z 

<p = —(tan 

2jt^ y-r d 


— tan 


-i 


y + r d 


) 


(3.15) 


where the origin of y coordinate is the center of the vortex pair and r d is half of the 
separation distance of the point vortex pair. After integation along the branch cut 
for the potential part (no contribution results from the circular contour of vortex), 
equation (3.12) becomes, using equation (3.13) for v 9 : 


n r ° ! • 1,-1 r * r » 

Di= ' , 2^ smh sr+'sJl \r.) ir 


(3.16) 


where r d is determined by using the conservation of centroid of the vorticity 


U = (R- R m ) jT‘ y^-dy = (R - R m )s, 


dy 


(3.17) 


By combining equations (3.1), (3. 4) and (3.10), equation (3.16 ) can be ex- 
pressed in terms of the tip aerodynamic parameters s t ,e t 


where e t is the tip aerodynamic efficiency 
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(3.19) 


e t 


C Lt 2 


it AR t Cpit 


Then, r 0) the radius of the turbulent vortex core, is expressed in closed form 
as below 

r„ = (*- W in -f (£) 2 7) < 3 ' 20 > 

Thus, the relation between the tip aerodynamic parameters; s t , e t , ( R — Rm) 
and r m (= r 0 ), and one of the vortex parameters, r 0 has been obtained in equation 
(3.20 ). Finally, r a will be determined by evaluating the integral in equation (3.20) 
as shown in the equation (3.47) in the following section. 

3.1.2 The turbulent viscous vortex core 

For large distances from the vortex center, the tip vortex behaves like a point 
vortex of constant strength r o . Near the center of the vortex , where the turbulent 
and viscous effects are significant, the vortex core can be described as a circular 
axisymmetric flow in space fixed cylindrical coordinate (r, 9 , z) as shown in Figure 
11. The assumptions are made that the tip vortex core structure is not affected 
by the inboard blade vortex sheet and the curvature of the vortex filament is 
neglected. 

With the further assumptions that the flow is incompressible with constant 
kinematic viscosity and the the axial gradient is smaller than the radial gradient, 
the flow in 6 direction in the vortex core can be expressed as follows : 


39 




u = u + u' 

v r = v r + v' r (3.24) 

Vg = Vg + v'g 


Substituting these velocities into equation (3.21 ) and using the continuity 
equation, we can obtain the equation for turbulent flow. Taking the mean and 
dropping for the mean values , equation (3.21) becomes 



(3.25) 
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where u is the kinematic viscosity v = /i/p. 

From the above equation, we can obtain the transport equation for circulation 
r in the core representing the angular momentum, rvg, as follows : 


d _ a_ ludv 

Dt r dr r dr 


~ 2 * + r £(«'^) 


(3.26) 


where 

T = 2t rrvg (3.27) 

The last two terms, including Reynolds stresses v' r v' e and u'v'g , represent the tur- 
bulent effects on the angular momentum in the core. The solution of the equation 
(3.26) can not be obtained without the knowledge of the two Reynolds stresses. 
However, an approximate solution which represents the main characteristics of the 
turbulent effect can be obtained if equation (3.26) is modified into an appropriate 
form such that the two terms including the Reynolds stresses are replaced by an 
eddy viscosity term as follows : 


DV_d_ 

Dt r dr 


. . iar 

^ + i) ~rTr 


(3.28) 


where e is the eddy viscosity in the vortex core which is assumed variable in the 
radial direction. 

The left hand side is the convective term and the right hand side is the viscous 
and turbulent diffusion term where the turbulent effect is much larger than the 
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viscous effect in high Reynolds Number flow. The eddy viscosity e, which has 
the dimensions of length times velocity, is modelled such that the length scale is 
assumed to be propotional to the radial distance r ; this is different from other free 
turbulent flows of constant length scale as pointed by Hoffmann and Joubert^. 
The velocity scale is assumed to be the radial gradient of the angular momentum. 
Thus, the eddy viscocity e may be written as 

e = k 2 r ~~~ (3.29) 


where A: is a constant which will be determined later in Section . 

The eddy viscosity tends to be zero toward the center of the core , where kinematic 
viscous is dominant, and also toward the outer turbulent vortex core boundary 
where the circulation assumes the constant T 0 of the potential flow. The eddy 
viscosity model is used throughout the turbulent viscous vortex core ; this form 
becomes the same as the eddy viscosity obtained by Roberts^! from the relations 
formulated by Hoffman and Joubert^l in a way analogous to the traditional 
mixing- length theory of Prandtl. 

When equation (3.29) is substituted into (3.28), it becomes 


ar 'ar ar a (, k 2 ai\iar 

m +U dl JrV 'Tr =T Tr ( ( " + r 2 hJ^-rTr 


(3.30) 


This equation must be solved for T subject to appropriate boundary conditions 
and using suitable approximate expressions for the convective velocities u and v r . 

In the axial direction, the velocity u varies due to viscous retardation and 
inviscid acceleration or deceleration such that the axial velocity u of the air just 
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behind the trailing edge of the moving blade tip would be the same velocity of 
the tip blade and it will return to ambient conditions (i.e., zero velocity) in a 
characteristic time T after passage of the blade. It is assumed that the vortex core 
persists up to time T (neglecting the initial time to roll up) and it will decay after. 
The average axial velocity in the plane normal to the centerline of the vortex core 
is approximated to vary proportional to a power of time, t , nondimensionalized by 
T , thus : 



(3.31) 


The radial velocity v r in the core is then determined from the continuity equation 
(3.23) after rewriting the axial variation in terms of the time such that x = Ut 


v r 


n r_ 

2 \ f) T 


(3.32) 


This inward radial veocity balances the outward turbulent viscous diffusion such 
that the vortex persists until the axial velocity becomes zero. The zero axial 
velocity after the time T results in no radial inflow with exponent n = 0 in (3.31) 
and (3.32). Then, there is only radial diffusion which causes a spreading of the 
vortex core and a decay in tangential velocity. 

With the above expressions for u and v T , equation (3.30 ) becomes, after rewrit- 
ing the axial variation in terms of time such that x = Ut, 



r dG 
\TJ T dr 


k 2 Y 0r d 

r dr 


dG. 1 dG\ 


(3.33) 
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where G(r,t ) is the nondimensional circulation in the core G = T/T 0 and c is a 
nondimensional parameter related to the inverse to the Reynolds Number, Re, 
based on the outer constant circulation T 0 in the vortex, : 

2jt v 2k 1 

C= ¥f^ = ¥Te 

Re is also related to the Reynolds Number Re t for which the characteristic length 
is based on the blade tip chord. By using the equation (3.4), 

Re = Re t • — (3.35) 

St 

where 

„ OR m 2(R-R m ) 

Rst — T7> 

v AR, 

(for a typical rotor blade c ~ 0.01 and Re ~ 10 5 ). 

The equation (3.33) can then be solved for nondimensional circulation G{r,t) 
in the core with the boundary conditions that the circulation is zero at the center 
of the core ; <j( 0, t) = 0, and the core merges with the potential solution, t.e., 
constant circulation r o , at the radial distance r 0 (t ) ; G(r 0 ,t) = 1. The initial 
condition for the decay period is given as the circulation profile at time T obtained 
in the persistence period. 

Approximated self similar solution 

The circulation profile T(t, r) is assumed similar. Then, the partial differential 
equation can be transformed into an ordinary differential equation by introducing 
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a self similar variable 


2 r 0 (i) r ! r °^ (r) 


(3.37) 


With this variable the equation can be written 


(c + zG ') — 
z 


+ 4/3 2 (7' = 0 


(3.38) 


where 



*Url{T) 

Ak 2 aT 0 T 


/ £ v2m+n-l 

(2 m + n) 


(3.39) 


To be similar, the parameter (3 should be independent of t and the exponents 
m and n should satisfy the relation 2m 4- n = 1 in the equation (3.39 ). In the 
persistence period, n = 1 , so m = 0 and r„(f) = r o (0) = r„(T). Similarly, in the 
decay period, n = 0 so m = 1/2 and r„(<) = r 0 (T)(t/T)i . 

The second order differential equation (3.38) , which is nonlinear, can be solved 
by an approximate method for small c with the boundary conditions 


G(0) =0 ;z = 0 (3.40) 

(7(1) =1 ;z = 1 (3.41) 


In the limiting case c = 0, the equation (3.38) is reduced to 


(G' 2 )' + 4(3 2 G' = 0 


(3.42) 
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The solution subject to the boundary conditions (3.40) and(3.41) is 

G = (1+ p 2 )z - p 2 z 2 (3.43) 

The approximate solution for c <!C 0 can be obtained after integration of equation 
(3.38) by seeking a solution in the form of a power series in small c. 

After integration, the equation (3.38) is reduced to 

— + G' 2 + 4 p 2 G = C x (3.44) 

z 

where Ci is constant to be determined later after applying the boundary condition 
at z = 1 in (3.41 ). The solution in the following form is tried by stretching the 
independent variable z and the dependent variable G for the nonlinear problem 
such that ; 

Q 

— = &>(?) + c + • • • (3.45) 


where 



(3.46) 


The independent variable z 2 instead of z is considered because of the behavior of 
the smooth variation of the even function charateristic of G across 2 = 0. The 
only valid value of the exponents p and q are p = 1 and q = 2. The solution G 
with the boundary condition at z = 0 in (3.40) is ( see Appendix B ) 
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(3.47) 


G=c(l'-l-ilny)-cM^j(C ll) -rini') 


where 

„ \/l + 4C 1 I» + 1 _ 2 2 

1 — . TJ = — 

2 ' ' c 2 

For 2 »c, near the outer boundary, this solution behaves as 

G » - /?V 


(3.48) 


(3.49) 


and for z <£. c, near the center of the core, 

z 2 

G « 2— (3.50) 

c 

The constant C\ is determined by appling the boundary condition G(l) = 1 in the 
equation (3.47), which becomes 

<?!«( 1 + P 2 ) 2 (3.51) 

Then, the solution for z c in (3.49) becomes the same as that obtained for 
c = 0 in equation (3.43). The arbitrary parameter (3 can be now determined by 
imposing the condition that the core smoothly merges with the potential solution 
G = G 0 ; G’(l) = 0, upon the solution in (3.47), giving 

/? = 1 + 0(c) (3.52) 
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The solution G(z ) obtained will be used to determine the turbulent core ra- 
dius r 0 in the equation (3.20) and the characteristics of the vortex core including 
the tangential velocity profile vg in the persistence and decay periods, and the 
persistence time T. 

3.1.3 The characteristics of the rotor vortex 

The main parameters that characterize the vortex velocity or circulation profiles 
are the radius of the turbulent core r 0 and the nondimensional quantity c (which 
is related to the inverse to the Reynolds Number), which are in turn related 
to the tip aerodynamic parameters ; the tip circulation shape function s ( , the 
tip aerodynamic efficiency e t , the maximum tip circulation T m and the distance 
from the tip to the position of the maximum circulation R — Rm- The additional 
parameters that characterize the vortex are the maximum velocity in the core V*, 
the radius of the maximum velocity r*, the circulation T* at the radius r* and the 
persistence time T. 

1) The turbulent radius r 0 is now determined from (3.20) and (3.47). 

In the persistence period , t < T, r 0 is constant with time 

r 0 (<) = r„(0) = r 0 (T) = [R - R m )s,/ sinh - h{c) j (3.53) 

In the decay period, t >T, r„ increases with time 

r 0 (<) = {R- Rm)st/ sinh - h(c)J ■ ^ * (3.54) 

where h(c ) is the nondimensional rotational energy in the core 
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(3.55) 


*M-r© 


2 


dr 

r 



G 2 — 
z 


The integral h(c ) is independent of the tip aerodynamic parameters but is generally 
a function of c because G obtained in (3.47) is a function of c such that the radius 
r 0 decreases with increasing the Reynolds Number. However, in the case of high 
Reynolds Number c« 1 the value can be approximatly obtained as by using 
the G in (3.49) and 0 = 1 in (3.51). Then, the vortex varies directly as ( R — Rm) 
which is a variable quantity depending on azimuthal blade position ( unlike the 
fixed wing) and it is also related to the shape of the tip circulation as s t ,e t . 

2) The tangential velocity of the vortex can be expressed as follows. 

In the persistence period, 


vg = 


2tt r 0 {T) 


V(z) 


(3.56) 


In the decay period, 


vg = 


2jt r a {T) 


™(?)‘ 


(3.57) 


where the velocity function V(z ) is defined as below by using G in (3.47) with 
/3 = 1 : 

V(z) = G(z)/z 

= c(Y - 1 - In Y/2) /z - z + c 2 (Y In Y) /Az , z< 1 (3.58) 

= 1 /z , z > 1 


where Y is given in (3.48) with C x = 4 and (3=1. 
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The velocity function can be approximated in a simple forms as shown in 


Figure 13 : 


V(z) = 2z/c , 
= 2-z , 
= 1 /* , 


Z « Z* 

2* < Z < 1 
2 > 1 


(3.59) 


where z* is the position of the maximum value of V(z). 

Near the center of the core, z « z‘, the viscous effect is dominant as in a the 
laminar sublayer. For z > z*, the turbulent effect is dominant up to z = 1 , and 
the vortex behaves as point vortex for z > 1. 

3) The location of maximum V is found by differentiation of V(z). For c<l, 



(3.60) 


and 

r* = r 0 (t)z * (3.61) 

The maximum velocity function V* is 

V-( 2 *) » 2 - (3.62) 

which shows that V* tends to be constant as c becomes small, that is, the Reynolds 
Number increases, whereas in laminar case V* increases as the Reynolds Number 
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increases. 


The maximum circulation function G* is 




ilal 

2 C 


(3.63) 


which shows, G* decreases as the Reynolds Number increases, whereas in the lam- 
inar case G* is constant. 

4) The persistence time T can be determined from equation (3.39) with /? = 1 


7rr o (0) 2 7rr o (0) 2 2 ARtS t 1 

4*T 0 = 4k 2 (R - R m )C Lt flRm 


(3.64) 


Thus, the number of revolution for which the tip vortex persist is: 


TCl _ 1 (R - R m )AR t sf 

2% 4A: 2 RmCu sinh 2 (4s 2 /e 2 — h(c)) 


(3.65) 


which is a function of the shape of the tip circulation and the Reynolds Number 
(through c). 


In the decay period, the rate of the diffusion can be obtained by differentiation 


of r Q (t ) with t such as ; 


dr 0 {t) _ r 0 {T ) i 


(3.66) 


The constant k of the eddy viscosity in (3.29) can be obtained from the expres- 
sion of T with the assumptions that it is universal for the different flows and is 
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independent of Reynolds Number. In the fixed wing case, the persistent distances 
( d = TU) are measured in many references'’!^], Then, the correlation between 
T of equation (3.64) for high Reynolds Number and the experimental data gives 
k = 0.06. 

So far, the rotor tip vortex has been defined and described in terms of the tip 
aerodynamic parameters. In the next section, the interaction between the rotor 

vortex and the following blade will be described especially when the maximum 

/ 

acoustic pressure is observed. Typically, the interaction is observed in l| revo- 
lutions of the vortex as shown in Figure 7. In that case, the vortex can be in 
the persistence or decay period, depending on the tip aerodynamic parameters as 
shown in equation (3.65). 
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3.2 Unsteady Loads During the Interaction 


With the tip vortex defined in the previous section, the unsteady field due to 
the interaction with the following blade is analyzed in the situation where the tip 
vortex is aligned parallel to the following blade span as shown in Figure 7, for 
which the maximum acoustic pressure is observed. The interaction is simplified 
such that the vortex is straight and stationary and the blade is moving with the 
uniform velocity Uf, within the finite region of the interaction as shown in Figure 
14. Thus, the interaction can be considered as two dimensional at each section of 
the blade span, neglecting the tip effect. The unsteady flow during the interaction 
will be described under the assumption of no distortion of the vortex path or of the 
vortex core. Even though the vortex is turbulent and viscous, it is assumed that the 
flow around the airfoil is a potential flow during the interaction. These assumptions 
mean that the vortex is forced, not free ; the vorticity field of the vortex core 
and the disturbed irrotational field, including the flow around the blade, and the 
acoustic pressure in the far field are coupled only through the surface boundary of 
the blade as explained more generally by KovcLsznay^l and Goldstein^. This is the 
same assumption as for traditional gust theory^. Thus the vortex behaves like a 
convecting wave. These concepts were visualized by Rogler^^' ,as shown in Figure 
3, in the problem of the interaction of a vortex array, representing a turbulent 
flow, with a semi-infinite plate. The stream function representing the vortical 
fields is not distorted but the stream function when combined with the disturbed 
irrotational field, satisfying the boundary condition of no flow through the solid 
boundary, is distorted to give a vortex on either side of the plate. This flow pattern 
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is in general agreement with the experiments of Yu^O) in the extreme case for 
which the center of the vortex core meets the blade. Within these assumptions, 
only the normal component of the vortex velocity on the blade surface is important, 
especially in the case of zero angle of attack between the convecting flow and blade. 
In general, the assumptions are valid when the normal component of the velocity is 
small compared to the convecting velocity and the rate of the diffusion in equation ' 
(3.66) in Section 3.1.3 is less than the convecting velocity ; t.e. when the interaction 
occurs in the diffusion period of the vortex. The flow near the blade is assumed 
incompressible in calculating the unsteady pressure on the blade surface in this 
section whereas the air is assumed compressible to obtain the acoustic pressure at 
the far field in the next section. 

3.2.1 Unsteady loading on the surface 

The governing equation for the incompressible small disturbance potential fuction, 

<f > , in the blade fixed coordinate is 

V 2 <f> = 0 (3.67) 

The linearized boundary condition for the thin airfoil that there is no flow across 
the surface becomes as below 

dd> 

— — -Wg{x - U b T) | x |< 6 2 = 0 (3.68) 

where w g is the normal component at the blade surface of the velocity of the vortex 
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convecting with the velocity U b and r is time. The combined velocity w g and the 
normal component of the disturbed velocity at the surface becomes zero. 


The perturbed pressure is given by the linearized unsteady Bernoulli equation 



'dj_ d£ 
8t 6 dx , 


(3.69) 


and the unsteady lift can be obtained by integration of the pressure difference 
between the upper and the lower surface along the chord 

L(r) = - / pu(x, r ) - p L (x, r)dx (3.70) 

J —b 

To evaluate the lift due ot the vortex interaction, it is useful to consider an 
oscillating blade in the uniform flow with the boundary condition that the flow at 
the surface moves with the same velocity of the blade : 


dj_ 

dz 


dzi 

dr 


+ U„ 


dz b 

dx 


x \< b z — 0 


= w b (x,r) 


(3.71) 


where 2 j(x, r) is the motion of the blade having zero thickness and w b (x, r) is 
a convenient shorthand for equation (3.71). For a simple harmonic motion of a 
blade, 
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A x)e tur 


z t (x,T) = 

Wt(z,T ) = Wb(x)e' UT (3.72) 


Using the conservation of circulation in the flow field (i.e. the shedding vorticies 
in the wake and the changing circulation around the blade) and assuming that 
the vorticies convect with the free stream velocity C/j , the pressure distribution 
on the blade in the simple harmonic case was obtained by Schwarz^! in the 
following form ( which can be obtained also by using the conditions of no pressure 
discontinuity across the wake, and at the trailing edge, i.e. Kutta condition) : 


-{Pu-Pl) 


2 p 0 






(3.73) 


where A p is the spatial part of —(pu ~Pl) and x* is the nondimension alized space 
xjb ,and k is the reduced frequency of oscillation defined as 



(3.74) 


From the result for the simple harmonic motion of the blade, the loading due 
to the gust can be obtained by comparing the boundary condition on the surface 
in (3.68) and (3.71). First, for a sinusoidal gust, the boundary condition becomes 
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W b = -Wg(x - U b T ) 


= -w g e 


tu(T-x/Ui,) 


(3.75) 


where w g is the magnitude and u is the frequency of the sinusoidal gust convecting 
with the velocity U b . Comparing (3.72) and (3.75), 

wi = — w^e~ ,kx ' (3.76) 

The pressure distribution on the blade due to the sinusoidal gust can then be 
obtained from (3.73) by using the above relation for wi. 

Next, for a sharp-edged gust (step gust), known as the Kiissner problem, the 
lift can be obtained by taking the Fourier transform of w g (x - U b T ) of the step 
gust^l, which is 


u> ff (w) = 



w g (x — U b T)e ,ljr dT 


—e~ ,k e~ ,k ** 

tu 


(3.77) 


where w 0 is the magnitude of the step gust. 

Then, 

wi = -Wg(uj) (3.78) 

and from (3.70) and (3.73) , the lift for the step gust becomes 
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L(t) = f \b f Ap(pUt,k,x*,-w g (u)dx t e VJT du (3.79) 

2ir J - co L J-i 


and may be expressed as follows 


L(s) = 2tt p 0 U b bw 0 ip(s) 


(3.80) 


where 


1 f a 


S(k)e ik l -V 


ik 


dk 


(3.81) 


which is the lift response function for the step gust, called the Kiissner functional 
and 


S(k) = C{k)[J 0 {k) - iJ^k)} + J x (k) (3.82) 


which is the lift response function due to the sinusoidal gust, called the Sears 
function^!’!® 1] and 


C(k) = 


H?(k)+xH^(k) 


(3.83) 


which is related to the lift due to the oscillating blade in the uniform flow, called 
the Theodorsen functional. 

J„ is the Bessel function of the first kind of order n and H ^ is the Hankel function 
of the second kind of order n. The variable s is the nondimensionalized time defined 
as 
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(3.84) 


For the vortex gust, the lift due to the interaction can then be obtained by 
using the Duhamel’s integral with the Kiissner function in the time domain. Thus, 


L(s) = 2npUb wJO)ip(s) + f — a) da 

Jo da 


(3.85) 


where a is the dimensionless time at the leading edge. 

Here, the normal velocity of the vortex gust on the surface, w g , becomes the 
tangential velocity of the vortex if the blade cuts through the center of the vortex 
core. After the transformation of the coordinate system between the vortex fixed 
coordinate and the blade fixed coordinate as shown in Figure 14, the velocity w g 


becomes 


_ y ( El 

2tt r 0 \ r a 


(3.86) 


and the gradient of the gust velocity is 


dwg(a) 

da 


_ r 0 t ' y, ( Xp 

2 nr 0 r 0 \ r 0 


(3.87) 


where V is the vortex velocity function of (3.58) in Section 3.1.3 and x 0 is the 
initial distance between the center of the vortex and the leading edge of the blade. 
As shown in Figure 17, the gradient of the gust velocity depends on the turbulent 
vortex size, r l} the laminar sub core radius, r* and the gradient at the center of 
the core 2/c. 
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The lift in (3.85) can be evaluated easily analytically or numerically if the 
Kiissner function, which needs an integration for the whole range of frequency as 
shown in equation (3.81), can be approximated in a simple analytical form. 

3.2.2 Behavior of the Kiissner Function 

The function is shown in Figure 15. It shows that for this step gust there is a 
time lag in the lift to have the asymtotic value. This time lag results from the 
wake contribution to the lift and the virtual mass effect due to the acceleration of 
the flow in addition to the effective change of the angle of attack due to the gust. 

The approximate form of the Kiissner function can be expressed in a series of 
exponential terms^l as follows 


rP(s) = 1 - 0.5e~°' 13 * - 0.5e“ 


(3.88) 


This is valid for the whole range of s in general. A more accurate form can be 
obtained by using the asymtotic expression of the function near the original, 
s < 2, (the time required for the front of the step gust passing the chord), as 
follows 


rp{s) 



s s 2 23 

j. — 5 I 

12 96 13440 J 


s< 2 


(3.89) 


Comparing (3.88) with (3.89) near the origin, there is no significant difference in 
unsteady loading as shown in Figure 15. But if this loading acts as the acoustic 
source, a significant difference occurs near the origin because it involves the gra- 
dient of the Kiissner function which might be called the acoustic response 
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function for a step gust. 


Then for s < 2, 

dtp(s) _ f2s /II 3 1 5 1 2 7 23 _ 2 \ 

ds -\ir\2s~ 2 12 + 2 96 S 2 13440 / 


(3.90) 


and for s > 2, 

d M s \ - o.5 • O.13e -013( * -01) + 0.5e -( * _01) 
ds 


(3.91) 


The gradient of the fuction has an integrable singularity at the origin and the 
expontial terms are modified to give a smooth function. This is shown in Figure 
16. 

The exact form of the gradient of the function is from equation (3.81) 

dJ aT = h II s <- k V l, '~" ik < 3 - 92 ) 

3.2.3 Evaluation of the lift 

The lift in (3.85) can be easily evaluated in the time domain with the approximated 
Kussner function in the previous section and the vortex velocity function V in 
Section 3.1.3. 


Analytic Evaluation 

First , with the simplest form of the Kussner function of (3.88) 


rp(s) = 1 - 0.5e -0 ' 13 * - 0.5e~' 
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and the approximated velocity function (3.59) in Section 3.1.3, the integral (3.85) 
can be calculated analytically. The velocity function V is transformed in the time 
domain with the relation used in (3.86), i.e., 


~ _ r _ x 0 - U\,t _ x Q - sb 
~ r a ~ r 0 ■ r a 

where x 0 is initial distance between the vortex center and the leading edge. 

The transformed approximate velocity function V'(s) is as shown below in each 
of five time periods. The velocity in the turbulent region and in the sublayer region 
on the vortex core are approximated by linear functions. 


(3.93) 


V.W = *5 

s 0 < S < Sx 

(3.94) 

V 2 (s) = -2 + ^ 

Si ^ s s? 

(3.95) 

v.M = r-ir 

S 2 ^ ^ S 3 

(3.96) 

V,(s) = 2+^ 

S 3 < S < S 4 

(3.97) 

r.(«) = 

s 4 < s 

(3.98) 


where 


Xq t 0 x 0 -r 

So — 0, Sx = , S 2 — 7 > s 3 


x 0 + r* x 0 + r 0 

L > S * ~ L 


(3.99) 


Because of the discontinuity in the gradient of the above approximated velocity 
function, the integration in (3.85) should be calculated for each period and it 
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should be added up to s to get the accumulated effect. Thus, the lift at s in 
s,_x < s < Si becomes : 


rp(s — a)dcx 

n — 1 

(3.100) 

The integration in each period is simple because of the constant gradient of the 
velocity except for the potential region which reqires the Exponential Integration 
functions i?i(a) and P,(a) as shown in Appendix C. The general form of the lift 
can be expressed as below 

^( S ) = 2xV^^ 5 j[P 1 W+P ^ We- I '■' 3 ’ + P3We-•] (3.101) 

where P,(s ), i = 1,2,3 are tabulated in Table 1. 

The lift becomes zero as s — * oo, which can be checked with the asymtotic 
value of the Exponential Integration function in P 2 (s ) and Pz{s) in s > s 4 . As 
shown in Figure 18, P,(a) diverges as s — > oo in the form : 


L(s) = 2 irtJbp 




rp(s — <j)d<j + 





a > 0 


(3.102) 





(3.103) 


whereas 1(a) converges. So, the expression for the lift becomes zero 
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e -p{'-¥) Ei jp(s - y)] -► 0 s — * oo (3.104) 

where p is positive. 

For the case of a discontinuity in the velocity function, the gradient becomes in- 
finite. In that case, the integration should be evaluated by the superposition of the 
regions having finite gradient with proper initial conditions at each discontinuous 
region. 

Numerical Evaluation 

A more accurate evaluation of the lift can be done numerically with the more 
refined Kiissner form of (3.89) for s < 2 and the complete form of the vortex 
velocity function V of (3.58) in Section 3.1.3. The integration in (3.85) should be 
calculated at each time s in the a domain. As shown in Figure 17, it is convenient 
to shift the coordinate a such that w g has a symmetric form. Then, the integration 
in (3.85) is equivalent to the following : 

Ms) ~ J k ^-{x)rp{s - x - j)dx (3.105) 

In discretized form, it becomes 

(3.106) 

(3.107) 


M S ) ~ 5Z - X n - T^ Xn 

n = 1 ax 0 


where 


x n = -^r + J2 A x i 

b ,= i 
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*">(_* i = *£» fl v 

dx 1 Si l " ) 


(3.108) 


To calculate the integral more effectively, the coordinate z is discretized nonuni- 
formly such that near the sublayer region of the core Az is small and in the poten- 
tial region Az is large. For the large intial distance (z 0 — * oo), the initial condition 
in (3.85) is neglected (tu fl (0) ~ ^ — » 0). In that case, numerical calculation is done 
for finite z 0 and the residual part (-oo < z < -z 0 ) is done analytically , where 
the velocity function (V = i) is that for the potential region. Then, the residual 
part of the lift is 


L R {s) ~ / x°° ““j I 1 “ 0.5e- 013( * +, - i,o/ ‘ ) - 0.5e- (l+ '- iro/6) ] dx 


(3.109) 


It becomes 


Ms) -77+O.5 — -77 — 0.13e°' 1S:ro/,4 i?i(0.13-j^) e-° 13, +0.5 -i- - e^E^) e" 

Xq/o Xq/u b XQ/b b 


(3.110) 


where 


r 00 g * 

E^k) = / — dx 

Jk x 

e~ k T 1 12 12 

~ —j— 1 H 1- — H 5- 

k L a cr a 3 


(3.111) 


(3.112) 
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3.3 Acoustic Pressure 


The acoustic pressure for a moving body in the presence of the vortex can be 
obtained by using the Ff owes- Williams and Hawkings equation in the medium at 
rest and by assuming the nonuniformity is confined to the sources ( i.e., constant 
acoustic propagating speed). The signal is analyzed in space-fixed coordinates and 
the sources are described in body-fixed coordinates. The equation has monopole, 
dipole and quadrupole sources corresponding to thickness, surface pressure fluctu- 
ations and unsteady flow fields. Here we consider only the dipole source because, 
assuming a thin airfoil and low mach number, we can neglect the monopole and 
quarupole terms. As explained by Curle^l f the dipole source in the equation 
represents the effect of the scattering due to the body for long wave lengths of the 
signal compared with the body size ; then the source can be related to the total 
force on the surface. Actually the observed acoustic signature during the blade - 
vortex interaction contains all frequencies but considering that most of the energy 
comes from the lower frequencies^' a chordwise compact source can give the 
basic property of the signal. In the case of a large but finite source region where 
the tip vortex and the following blade interact , a spanwise noncompact source 
is more appropriate. The acoustic pressure is analyzed in the simplified situation 
such that a finite unsteady dipole source, whose strength is related to the unsteady 
load, moves linearly. 

The inhomogenous wave equation with a moving dipole source is 

(£-* v, )('-* ) --£H')£) (3 - n3 > 
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where a D is the constant speed of sound in the medium at rest, F tJ is the dipole 
strength defined in (2.11), and / is a function of the body surface. 

The solution of the acoustic pressure p in 3-D space is 


p(x, t) 


—— [ 
4?r dx x Js 


F *j n j 

R | 1 -M - f 


dS(f}) 


(3.114) 


where p is related to p with the isentropic relation p = a? 0 {p — p 0 ). 

The relative distance between the observer fixed and the source in linear motion 
with the velocity U is 


R =| x — i Ur e — rj 


(3.115) 


where x„ is the initial observer location in the space-fixed coordinates and rj is the 
source positions in the blade-fixed coordinate as shown in Figure 19. 

x — R\i -J- R?j -f- R^k (3.116) 

fj = f?i» + Vij + r) 3 k (3.117) 

The source is evaluated in emission time r e defined as 


, i2(r e ,£,^) 

Tc = * ~ 

a 0 


(3.118) 


The above time relation is related to the Doppler shift due to the relative motion 
between the source ( rj ) and the observer (i) and 


1 

1-M R/R | 


(3.119) 
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is the Doppler factor in amplitude due to the moving source. 

3.3.1 Far field approximation 

The solution can be expressed in terms of a time derivative instead of a space 
derivative by using the relation of (3.118). For the far field, it can be approximated 
as in equation (2.57) 


p(i ' t) = -4kL 


V A d "(9) ( 3 - 120 ) 


where /,• is the force on the blade. 

Neglecting the time variation of R along the source region and the drag force 
effect, the acoustic pressure can be written as follows for a chordwise compact, 
but spanwise noncompact, source 


p(x,t) = 


4ira 0 £* 2 | 1 - M • | 


r 

-ltr.ni i2 J- 


+h 

h 


Tr L(T ’ m) 


dr) 3 


(3.121) 


J r e 


where R* is the distance from the observer to the center of the source and / 3 is 
half of the source distance in the spanwise direction rj 3 as shown in Figure 19. For 
a rectilinear moving source, the emission time r e is obtained as shown below : 


l _ SjM. \j(t- - (1 - M 2 ) { t 2 - 

1 -M 2 T 1 -M 2 


(3.122) 


The + sign is for the source moving toward the observer. 

The acoustic pressure is composed of two parts; a directivity function for a 
point source defined as 
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(3.123) 


R* | 1 - M • | 2 

and an integral along the spanwise source direction for a noncompact source. 

The integrand of the time variation of the lift in the integral can be obtained 
as below by using the Duhamel integration of the gradient of the Kiissner function 
defined in (3.90) and(3.91) in Section 3.2.2, ( which represents the acoustic 
response function for the step gust) 


d L 

— (s,y 3 ) = 2irpU b b 


. drp(s) 


/; 


dwg(a) rp(s — a ) 


da 


ds 


dal 


(3.124) 


The acoustic pressure also depends on the gradient of the gust ^ which is a 
function of r Q and c as shown in Figure 17. The schematic diagram for the 
Duhamel integration is shown in Figure 20, where the contribution of the gradient 
at the center of the core is large because of the delta-function-like acoustic response 
function The contribution to the integral along the span of equation (3.121) 
comes from two factors: one from the lift variation along the span, which is zero 
in our 2-D approximation in the calculation of the unsteady loads , and the other 
from the emission time variation along the span as seen in equation (3.118); 

3.3.2 An alternative form for a spanwise noncompact 
source 


The acoustic pressure for a noncompact source expressed in the form of (3.121) 
can be interpreted such that the acoustic pressure observed at x and t is the 
accumulated pressure of the multiple compact point sources emitted at the same 
time r e but having different propagation times R(r e ,x,Tf)/a 0 depending on the 


70 


position of the sources fj as in (3.118). And it has an integral for which the 
variable r) 3 is different from the variable of the integrand of the gradient of the 
unsteady lift expressed in term of r e The difference in the variables makes the 
integration sometimes difficult to evaluate. An alternative form of the acoustic 
pressure!*^ can be considered which has the same variables for the integrand and 
the integral as follows : 

For fixed x and t, drj 3 in equation (3.121) can be written 


drj 3 — 



(3.125) 


Then, equation (3.121) can be written, after changing the varable rj 3 to r e , as 


p(Z, t) 


1 r e+ <L l(t V ___R L ___djT 1 
4na 0 Jr~ dr K eJ R? | j _ M ■ § [ 2 dr e 


(3.126) 


where r e + and r e are r e at rj 3 = l z and —hi respectively in equation (3.122). 

The quantity ^ can be obtained by differentiating the time relation (3.118) with 
r e and using (3.115), i.e., 


1 + 


J _dR 
a 0 dr e 


= 0 


(3.127) 


where 


dR 

dr t 


R f d n TT \ 


(3.128) 


Then, 
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(3.129) 


* |i = 1 * £ 

R dr. R 


and becomes 


dy 3 R( 1 - f • M) 

ar, R3-V3 


(3.130) 


Thus, the acoustic pressure p(z, t ) can be written as 


^ 1 ( Tt d r, X ^ 

p(x,t) = -- / — L(r e ) ^ — a 1 

4;ra 0 J T - dr j? | ! - M • § | (fl 3 - > 73 ) 


(3.131) 


where 


R ■ M = — Ut,r e 

R 3 -V> = v/* 2 - W - *V«) 2 - % 

t-T. 


(3.132) 

(3.133) 

(3.134) 


Neglecting the time variation in the point source directivity part as in equation 
(3.121) for the far field , it becomes 


1 ^2 T' d r ( a 

^ y) 


' f » -fM) 


R* - vs 


dT e (3.135) 


With nondimensional variables, 


1 R> ft d , (l-hP^js-t + R.M 2 

- — / — — L(sJ ===== — — as 

Ta o R* 2 I 1 - M • I 2 My (< - s) 2 - M 2 ((Ri - s) 2 + R?) 


. 4jtc 


M\/{i - s) 2 - M 2 ((fl x - s) 2 + J2/) 

(3.136) 
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where R t and Ri are nondimensionalized by the semichord 6, and t is nondimen- 
sionalized by b/Uf, and 


7_e ta ,/(? - 2£e) 2 - ( 1 - m=) {? - 

5* = =F 


1 -M 2 


1 - M 2 


(3.137) 


The acoustic pressure expressed in (3.131) or (3.135) can be interpreted such that 
the the acoustic pressure observed at x and t is the accumulated pressure emitted 
at different times r e depending on the source position y 3 but having the same 
propagation time R(r e , x, fj)/a 0 . 

The forms derived here are useful in the evaluation of the integral because the 
integrand is a direct function of the variable s. It will be used in the calculation 
of the acoustic pressure especially when the observer is far off the center of the 
source. In that position, the observer can easily detect the effect of the spanwise 
noncompactness of the source having relatively large aspect ratio. 

3.3.3 Interference of the spanwise non-compact source 

The compactness and the noncompactness depend on the characteristic time of 
the acoustic signal and the source dimensions. In most cases, there is no clear 
relationship between the charateristic time of the acoustic signal, which is not a 
simple harmonic signal, and the source dimensions whether the source is compact 
or noncompact. To visualize this relation, it is useful to consider a point source 
signal and to interpret the signal due to the noncompact source as the result of 
the interference of multiple point source signals. 
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Consider a point source signal, first, which can be represented as dL(s)fds in 
the emission time domain s as shown in Figure 21 neglecting the point source 
directivity function in equation (3.123). In the emission time domain, the lift L(s ) 
can be the one measured or computed for each span location. The acoustic signal 
really observed is the one in the reception time domain ? as shown in Figure 21(b) 
which can be obtained after transformation from the emission time domain with 
the time relation t = r e +R(r e , rji)/a 0 . The amplitude is the same for corresponding 
t and r but the width d t reduced for the source moving toward the observer. Then 
the noncompact source signal can be obtained by summing the point source signals 
in the reception time domain t with the delayed time A t due to the different source 
positions y for each signal as shown in Figure 21(c). This delayed time related to 
the source dimension causes interference, which contributes to the amplitude in 
a destructive way and alsp changes the directivity pattern from the point source 
directivity. 

The interference is closely related to the ratio of At and d t , where the width 
of the signal, d t , depends on the strength and the strucure of the vortex and the 
maximum delayed time between the ends of the source region, A t max , depends on 
the observer position, as shown in Figure 22. If the observer is far from the plane 
of symmetry of the source, greater interference is expected, but if the observer 
is in the plane of symmetry of the source, there is small interference where the 
maximum acoustic pressure is observed ; then the source can be approximated as 
a point source to evaluate the maximum acoustic pressure in this simple geometry. 

The effective ratio of At and d t can be estimated from a simplified point source 
signal and its spectrum as shown in Figure 23. Assuming the signal is periodic, it 
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can be expanded into a Fourier series having a basic period, T 0 and a fundamental 
frequency, u 0 = 2ir/T 0 . Then, the signal /(£) can be written as 

/(£)= £ F(nu o)e• Wo, (3.138) 

n=— oo 

where 


1 f T °/ 2 

n™ o) = w , /(()£-"""'* 

T 0 J-To/2 


(3.139) 


For the simplified signal of a generalized sawtooth function, the Fourier coefficient 
■F(nu/ 0 ) is 


F(nu 0 ) = 


1 - e-» 2rn * - t2irnp, 
(2nnp t ) 2 


■ Ap t 


(3.140) 


where A is the peak amplitude and p t is the solidity defined as d t /T 0 . As shown 
in Figure 23, the amplitude of the nth harmonic | F(nu 0 ) \ decreases toward the 
higher harmonic as ( np t )~ 2 . 

Assuming that the spectrum begins to decay significantly at np t — N, the 
corresponding harmonic is 

N 

n e = — (3.141) 

Pi 

which becomes higher for smaller p t , that is, narrower the width d t . The corre- 
sponding frequency is 




n c u>o 

N 

— U>o 
Pt 


(3.142) 
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Thus, the corresponding period is 


T = 

J-ne — 


2?r 
w ne 

—El 

N u 0 
dt_ 

N 


(3.143) 


This period T ne can be the characteristic time of the acoustic signal and the con- 
tributions of the waves for. which periods are less than T nc is small. From Figure 
23, it is seen that N = 1. Thus from equation (3.143), T nc = d t . 

Thus, if d t is much greater than the time due to the delay in the source, 
At, there will be no interference between the signals as for a point source. Or 
interference will occur approximately when A t > d t /A. 

3.3.4 Evaluation of the Acoustic Pressure 

For the noncompact source or the point source, the acoustic pressure depends on 
the lift variation with time dL/ds due to the vortex interaction. 


Point source 

For the point source, the acoustic pressure can be described as shown below from 
equation (3.121) 


1 M 2f 3 dL(s) 
p = -- — — — D 


4wa 0 R* b 


ds 


(3.144) 


The directivity function D expressed in terms of directional angle is 
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where 


D = 


sin 9 

1 — M cos# cos $ | 2 


(3.145) 


Ri = R*cos9cos<f> 

R 2 = R*sm9 (3.146) 

R 3 = .R* cos# sin 0 

This directivity function evaluates the effect of the observer position located in 
the far field with the same radius from the initial source position rather than from 
the moving source position as shown in Figure 33. 

From the general expression for lift L(s) in equation (3.101) in Section 3.2.3, 
the lift variation can be obtained as 

~L(s) = ^{P,( S ) + / , 2 (s)e-°' ,3 ‘ + P 3 (s)e-} 

= Q,(s) + Q 2 (s)e-° AS - + Q,[s)e- (3.147) 

where Qi(s),i = 1,2,3 are tabulated in Table 2. 

From this expression, the maximum acoustic pressure can be obtained by dif- 
ferentiating the pressure with time. Differentiating -§^L(s) between s 2 and s 3 , 
the pressure increases monotonically and differentiating it between s 3 and s 4 , the 
pressure decreases monotonically. So the maximum pressure occures at s = s 3 and 
the maximum pressure in nondimensional form is 
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(3.148) 


+ zo^iy*) (e-*^ -i) 

- £>.5 k i -k i 'fE,(k i rJb)e-'“T 

1=1 


where 


r , N ,( 1 M2l 3 pU b bT 0 

/w = - p/ Uft ~tt d 


(3.149) 


and k x = 0.13 and A: 2 = 1. 

For the extreme case c = 0, (after series expansion of terms including for 
small c) 


= -^-Eo.sfe + ( 3 - 15 °) 


which gives a finite value of the maximum acoustic pressure. 

As mentioned before, a more accurate results can be obtained by using the 
acoustic response function of (3.90) in Section 3.2.3 for s < 2. In this case, we 
need numerical integration as in the case for the evaluation of the lift of (3.106) 
in Section 3.2.3 except near s = x 0 /6 which has a singularity in acoustic function. 
This can be done by analytic integration for s < e with the linearization of the 
gradient of the vortex velocity ^ for small e and numerical integration for s > e 
as below 


dL{s ) A dw g dip x 0 

-JT = I ‘ + £ ' aJ (s - - 7 )Ai " 


(3.151) 
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where 


Zi 

II 

1 

Hi? 


Xn 



X„ 

= -He a*, 

(3.152) 

e 

0 i=i 

< Axn 



The analytic part I ( is 


T _ /■— Tp- V2/1J_ 

' \2^y 


1 5 1 7 23 

_ /y _| y 3 / 2 

8 V ^ 2 96* 213440 



dx 


(3.153) 


where a and b' are the coefficient of the linearized velocity function such that 


w g « ax + b' 


and 

u>„(s - f + e) - 10,(5 - f) 

CL — " 

e 

6 ' = a{-s + ^-) +w 9 (s -j) 

y = 3 x -j— (3.154) 

o 

The residual problem for x 0 — ♦ oo can be neglected as in the case of the evaluation 
of the unsteady lift in (3.110) in Section 3.2.3. 

With this more accurate acoustic function it turns out that for the extreme 
case of c = 0 the acoustic pressure has peak value of infinity. This is different 
result from the finite value of (3.150) which is obtained by using the approximate 
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acoustic function. The infinite value is due to the discontinuity in the vortex 
velocity w g for the case c = 0, in which case the velocity can be obtained by 
superposition of the simple velocity function including the sharp edged gust. The 
corresponding acoustic pressure to the sharp edged gust is the acoustic function 
defined as (3.90) in Section 3.2.2 which has the infinity at the origin. Thus, the 
acoustic pressure for this case c = 0 has the infinite peak value. 

The signal evaluated is in the emission time domain s and its numerical values 
are obtained at nonuniformly divided intervals of s as in the case for the evaluation 
of the unsteady lift to calculate efficiently. The signal in the reception time domain 
t is then obtained after the transformation of the numerical value of the pressure 
by using the time relation : 

' i= 

at 


Noncompact source 

From the results of the point source, the pressure for the noncompact source can 
be obtained from equation (3.121) or (3.136) depending on the observer position. 
If the observer is near the plane of symmetry of the source, equation (3.121) is 
more efficient in calculation because the results converge after the summation of 
relatively few delayed point source signals. The process of the calculation is the 
same as shown in Figure 21. The practical problem comes from the fact that each 
point source signal in the reception domain 1 has numerical values at nonuniform 
time intervals of t ( due to the nonuniform intervals in the emission time domain s 
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and its transformation to the reception time domain through its own time relation 
depending on the source position ). This causes the difficulty in numerical sum- 
mation of the point source signals. The problem can be solved by reevaluation 
of the numerical values of the point source at uniformly distributed values of ? 
close enough to pick the peak value of the signal. This requires the interpolation 
of the values nonuniformly distributed . Once it is done, the summation of the 
delayed point source signals is easy for the case that the results converge quickly 
with relatively few point source signals. 

For the case which requires many source. signals to converge, the other form 
of (3.136) is more efficient in the evaluation of the acoustic pressure as in the 
case where the observer is far from the plane of the symmetry. In this case, the 
accumulated acoustic pressure at each value of t , coming from different source po- 
sitions and emission times, is evaluated by integration of the integrand from sj to 
s+ numerically in the emission time domain s, so it does not require interpolation 
to sum the signals. This method has its own problem especially when there is a 
small interference effect between the sources ( In order to include the peak value 
of the acoustic pressure the pressure should be evaluated at each small range of t, 
which 'requires a lot of integration ). 

In both methods, to get the directivity effect of the peak acoustic pressure 
due to the interference of the noncompact source in addition to the point source 
directivity function D of the equation (3.145) , the geometric relations between 
the distance R and the directional angles in (3.146) are used in the time relation 
of (3.155) for the time transformation and in the integrand of (3.136). 
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3.3.5 Moving observer and moving source 

So far only for the moving source, which cause the Doppler shift and Doppler 
factor, has been discussed. If the observer is also moving, the time relation related 
to the Doppler shift in (3.118) is changed, that is, the relative distance R between 
the observer and the source is changed according to 


where 



R{T e ,t,x,y) 

do 


(3.156) 


R — | (Ro -f- V t — U T e )i + R 2 j + R$k (3.157) 


and V is the velocity of the observer moving at the same direction of the source. 
So the time relation of (3.155) is changed such as 


t = r e ( 1 


M 0 M) + M 0 — 
a 0 


+ 


Te( 1 - M 0 M) + M 0 M 2 r e 2 - + 52+5i 1 


\ 


1 -Ml 


1-M 2 


(3.158) 


and s+ and s~ in the equation (3.136) are changed after substitution of Rq + Vt 
instead of R x in equation (3.137). 

The above time relation is reduced to R 2 = R 3 — 0 to 


(l-M)r e + % 

(1 -M 0 ) 


(3.159) 
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, which becomes "the same as (2.34) with Ci = — 

The acoustic pressure for a point source or a noncompact source can then 
be evaluated for a moving source and a moving observer by following the same 
procedure of the previous subsection with the time relation (3.158). In this case, 
the directivity function can determine the effect of the observer position which is 
moving at the same velocity of the source, that is, located at the same radius from 
the moving source rather than from the initial position of the source. 
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3.4 Noise Reduction and Performance 

From the previous section, it is known that the acoustic pressure at the far field 
depends on the delta-like acoustic response function and the gradient of the vortex 
velocity at the surface foF given source and observer motion. The gradient is 
directly related to the vortex size r 0 , which can be expressed in terms of the blade 
tip aerodynamic parameters, and is also related to the Reynolds number Re via the 
nondimensional parameter c. In general, the amplitude of the acoustic pressure 
is reduced for larger vortex size and lower Reynolds Number. However, a lower 
Reynolds Number corresponds to a lower maximum circulation along the blade 
span r m , which is one of the primary parameters in blade performance. So, with 
fixed Reynolds Number, the acoustic pressure may be reduced by increasing the 
vortex size. 

As illustrated in Figure 24, several methods can be considered! to increase 
the vortex size.- The first group of methods is the modification of the tip shape, 
the second group is the attachment of the vortex dissipator to the tip, and the 
third group is the injection of air from the tip. The injection of air at the tip 
along the center of the vortex accelerates the axial flow near the vortex core. At 
the same time, this also reduces the inward flow at the beginning of the formation 
of the vortex, which reduces the persistence time T. Then there is only outward 
diffusion near the core and when the vortex encounters the following blade after 
time i, the vortex size increases from the initial size r 0 (t = 0) to r o (0)\fT/t due 
to the diffusion as expressed in equation (3.54) in Section 3.1.3. Practically, the 
alignment of the direction of the injection with the center line of the vortex core 
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is not easy and this also require additional power to inject the air. The vortex 
dissipater considered in the second group also promotes the earlier diffusion of the 
vortex core. However, this increases the form drag at the tip. And the two kinds 
of the methods explained requires the additional analysis to describe the detail 
relationships between the parameters related to the amount of the air injection or 
the geometry of the dissipator and the vortex core. So the first group of methods 
( increasing the vortex size by modification of the tip shape ) is considered now 
because the relationship between the initial vortex size and the tip aerodynamic 
parameters is already defined in equation (3.53) and (3.55) in Section 3.1.3 as 

r o (0) = (R- R m )s t / sinh ^ j (3.160) 

As shown in the above equation, the vortex size increases with the distance between 
the position of the maximum circulation along the span and the blade tip, t.e., 
( R — R m ). However, the larger distance a loss of lift considering that the lift 
is weighted toward the tip due to the rotation of the blade. So the analysis is 
restricted for the given distance ( R — Rm) in addition to the given maximum 
circulation T m , which is mentioned before. Then the remaining parameters are 
the tip loading factor s t and the tip aerodynamic efficiency e ( , which are directly 
related to the circulation shape near the tip. Before further analysis, it is assumed 
that the tip circulation shape normalized by the maximum circulation Y m and the 
distance (R — R m ) are almost constant during one blade revolution. For forward 
flight of a helicopter with low advance ratio and tip shape fixed, there is a resonable 
approximation even though in practice and (R — Rm) depend on the azimuth 
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angle ip. 


3.4.1 Modification of the tip circulation shape - 

To make the anaysis simple, consider the effect of the modification of the tip cir- 
culation shape rather than the tip blade shape itself. And consider the normalized 
tip circulation shape expressed in terms of the angle 9 rather than y ( for which 
the origin is the position of the maximum circulation) as follows 


where 


and 



y = Y/(R - Rm) 
= cos 9 


o <y < i 

- > 9 > 0 
2 ~ “ 


(3.161) 

(3.162) 


(3.163) 


(3.164) 

(3.165) 


The general tip circulation shape should satisfy the following boundary condi- 
tions 


G(0) = 0 
*<¥> = 
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0 


(3.166) 

(3.167) 


The circulation shape considered here, which satisfies the above conditions, is the 
form of the power of the sine function : 


G(0) = (sin 0) p 


(3.168) 


This function represents the elliptic loading if p = 1 as shown in Figure 37(a) and 
the rectangular loading if p = 0. 

The loading factor s t is then obtained as follows 

s t = [ G{y)dy 

Jo 

= I 1 (sm 9) p+1 d9 (3.169) 

Jo 

If p is an odd integer 


1 ■ 3 • 5 • • • p 7T 
2-4-6 — (p H- 1) ’ 2 


(3.170) 


and if p is an even integer 


2 • 4 • 6 • • • p 
1 -3-5 •• (p+1) 


(3.171) 


Ifp is natural number rather than integer 


_ r(f + i)r(|) 
‘ zr(f + |) 


(3.172) 


where T(i) is the Gamma function having property 
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(3.173) 


r(z + 1) = ir(z) 

To calculate the aerodynamic efficiency e t , it is convenient to expand the cir- 
culation function (3.168) into Fourier sine series : 


G = (sin 9) p = Y A m sin m9 


(3.174) 


where m should be odd integer to satisfy the conditions of (3.166) and (3.167). 
The coefficients A m is obtained as follows 




sin kffdd 


(3.175) 


where 


/ * sin m6 sin k9d9 = 0 , m ^ k 

Jo 


7T 


= 4 < m = k 


(3.176) 

(3.177) 


If p is integer and odd number, the coefficient A m is obtained from (3.175) and 
(3.174) as follows 


(-1)^ 


2 p 


-l 


sinp0 — ( v sin(p — 2)9 + • • •( — 1)** 


! ( 40 sin * 


(3.178) 


(:)* 


where I is the binomial coefficients defined by 


(!) 


p(p-l)(p — 2)--(p — g + 1) 


(3.179) 
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For example, if p = 1, Aj = 1 and if p = 3, Ai = | and A z = — i. 

If p = 0; the coefficient A m can be directly obtained from equations (3.175) 
and (3.174) as shown below 

A m = — — (3.180) 

7T m 


If p is integer and even number or natural number, the coefficient A m is 


Am — 


4 r § 

— / (sin 9) p sin mddB 
n Jo 

TC JO 

A + II + 1$ + 


\P +1 


(2cos0) 


m — 1 


m i 2 j (2cos0) m ~ 3 + ^ m 2 3 j(2cos0) m - 5 - 

(3.181) 


where 


/, = -(-I)*" 1 ( m * ) [ 2 (sm0) p+l {2cos9) m -^-^ 
IT \ 1 — 1 ) Jo 

r ifiy-l f m “ * ^ om— 2>'+l r ( 2 + jj£(f - *+ *) 

7T l ' v 1 / 2r(^ -» + 2) 


(3.182) 


From the coefficients A m obtained, the aerodynamic efficiency e t is easily ob- 
tained from the well known expressions'll as below 


and 



(3.183) 


s< 



(3.184) 
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Typical values are e t = 1 and s t = f for elliptic loading and e t = 0 and s t = 1 
for rectangular loading. 

Thus the radius of the vortex size r 0 can be obtained from the equation (3.160) 
by using the above expressions of e t and s t for the normalized circulation shape 
(3.168) as below : 


r o(0) 

R — R-m 


— Ai/ sinh 
4 




(3.185) 


The nondimensionalized vortex size r 0 /(R — R m ) is 0.3488 for elliptic loading and 
0 for the extreme case of rectangular loading ( in which case the tip vortex is point 
vortex and the induced drag, which is related to the rotational energy of the tip 
vortex, is infinite). 

From the analysis done so far, the acoustic pressure due to the blade vortex 
interaction can be calculated for the tip circulation shape considered here, and 
the relation between the noise reduction and the corresponding tip aerodynamic 
efficiency, which is one of the parameters of the performance, can be obtained for 
the various circulation shapes, that is, for various powers of the sine function of 
equation (3.168). 
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Chapter 4 

Results Discussion 


From the previous analysis, we find that the unsteady lift and the acoustic pressure 
are closely related to the slope of the vortex velocity, which is a function of the 
vortex turbulent core size r Q and the nondimentional parameter c ( inverse to 
Reynolds Number Re). 

The vortex velocity profile has been computed and compared with the avail- 
able experimental data. Unsteady lift and acoustic pressure are calclulated and 
variation of maximum acoustic pressure is found in terms of the parameters r x and 
c. The noise/performance trade-off for blade tip loading shape is also discussed. 

Vortex Velocity Profile 

The general vortex velocity profile along the radial distance r from the center of 
the vortex core is shown in Figure 13, based on equation (3.58). The distance r 
and the velocity v& are nondimensionalized by the turbulent core radius r 0 and by' 
the potential velocity at r 0 , that is, r o / 27 rr OJ respectively. 

The specific vortex structure is calculated and compared with the experimental 
data obtained by Tungl^l and Ballard^] for a blade having AR = 13.7 and R = 
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1.05m in the hovering condition. The experiment results are shown in Figure 25 
and Figure 26, respectively. In Figure 25, Tung’s data measured by hot-wire show 
the fluctuating velocity after the peak value, indicating the turbulent behavior of 
the flow, and shows a smooth variation of the velocity in the potential region and 
near the center of the laminar core. The circulation in the potential region seems 
to be constant and its value is almost the same as the maximum circulation on the 
blade. In Figure 26 , Ballard’s data measured by a laser velocimeter are shown 
with the circulation along the blade span. 

First, Tung’s data is compared with the present theory and it is shown in 
Figure 27. Both the velocity and circulation profiles of the vortex show good 
agreements with experiment in general, even though the maximum velocity and 
its position are underpredicted. In the experimental data, the circulation of the 
vortex r o is almost 90% of the maximum circulation at the blade The velocity 
profile predicted with the present analysis is also compared with that predicted 
with the Betz and the Spreiter and Sacks analysis for the persistence region in 
Figure 2. In the case of Betz theory ( modified by Donaldson^), it shows 
generally good agreement outside of the core. Near the center of the core the 
present analysis is more realistically represents the turbulent viscous core whereas 
the Betz vortex has infinite velocity, especially for the elliptic tip loading. In the 
case of the Spreiter and Sacks theory, the maximum velocity is finite but the its 
position is overpredicted because of the solid rotation model of the core. 

It is assumed that the vortex at ip = 65.4° is in the persistence period based 
on the estimation from equation (3.65) that the vortex presists until rp = 150° for 
these conditions. And the distance from the location of the maximum circulation 


92 


to the tip, ( R — Rm), is taken as 4.4 % of the blade span length, R. The value is 
obtained from the comparision between the measured circulation near the center 
in Figure 25 ( which corresponds to G{r/R) = 130 • r/R ) and the simplified 
expression of the circulation ( G[z ) = V(z)z = 2 z - z 2 where z = r/r 0 and 
V(z) = 2 — z ) from equation (3.59). From the two expressions for the circulation 
, r a becomes 0.0154i2. Assuming the tip circulation is elliptic, r„ = 0.349(72 — i?^ 
from equation (3.53). Then (R — R m ) becomes 0.04472 by comparing the two 
expressions for r 0 . 

In terms of the blade geometry, the radius of the turbulent core r„ = 0.4226 
( where 6 is the half of the chord length) and the diameter of the core is 3.5 
times the maximum thickness of the NACA 0012 airfoil. The maximum velocity 
vg is almost 40% of the rotating tip speed Q72. Reynolds Number based on the 
circulation in the potential region, Re = T 0 /u, is 2.0xl0 5 . 

Next, Ballard’s data are compared with the simplified equation of the vortex 
in equation (3.59). The vortex just behind the trailing edge in Figure 26 seems 
to be not completely rolled up and the vortex at ip = 180° is assumed in the 
persistence period ( the vortex persists untill almost two blade revolutions from 
equation (3.65), which is longer than the Tung’s vortex because the maximum 
circulation is lower and the distance (R — R m ) is larger). Taking (R — R m ) = 0.1R 
from the Figure 26, r 0 = 0.035R for R c = 1.0x10 s . Then, the position of the 
maximum velocity r* = 0.18r o = 0.0071? from equation (3.60), which seems to be 
smaller than that obtained from the data. However, the trend in equation (3.59), 
that the radius r* becomes larger for a larger r 0 and a lower maximum circulation 
r o ( ».e., higher c ), is as expected comparing the data obtained by Tung with 
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the data by Ballard. The gradient of the velocity, nondimensionalized by the tip 
speed, in the turbulent region seems to be 3.75 in Figure 26 and is 3 from equation 
(3.59), and the maximum velocity is around 20% of the tip rotational speed both 
in the data and the simplified expression V(z) = 2 in equation (3.59) ; this is 
considered to be good agreement. 

From the comparisons with the experiment data for the two cases, the present 
theory seems to give a reasonable estimate on the velocity of the tip vortex. It 
should be also mentioned that the maximum velocity of the vortex is around 
20 ~ 40% of the rotating tip speed in these cases. But if the position of maximum 
circulation ( R — R m ) is shifted toward the root, the maximum velocity of the vortex 
can be reduced substantially. 

Unsteady Lift and Acoustic Signal 

Typical nondimentional unsteady lift and acoustic signal for the vortex gust are 
plotted in Figure 28. It is observed that the lift calculated from equation (3.101) 
is similar to the vortex velocity profile because the negative velocity of the vortex 
contributes to the change of the effective angle of attack at the blade surface, 
which results in a negative lift at the beginning of the interaction. It is shown 
that the vortex velocity is symmetric at the surface but the lift is unsymmetric 
due to the behavior of the Kiissner function as shown in Figure 15 ( It introduces 
a time lag near the origin of the function due to the shedding vortices and the 
apparent mass effect ). 

The acoustic signal shown in the Figure 28(b) is for the part of the lift variation 
in equation (3.147), which is unaffected by the location of the observer. The signal 
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is calculated for a point source in the emission time domain s and it is similar to 
the slope of the vortex velocity profile because the acoustic function in equation 
(3.90) acts like a delta function which assumes the value of the slope during the 
integration in equation (3.124). So the amplitude of the peak acoustic pressure is 
related to the slope of the vortex at the center of the laminar sub core, 2/c and the 
width of the signal is related to the radius of the laminar core size, r*. The signal 
is also distorted due to the time lag in the acoustic function as in the Kussner 
function. From the acoustic point of view, the unsteady lift acts as a strength 
of the source and its slope is directly related to the amplitude and width of the 
acoustic signal through the vortex structure. 

Reynolds Number and Vortex Size Effect 

It is expected that both the peak unsteady lift and acoustic pressure decrease with 
increasing turbulent vortex core size r 0 and smaller Reynolds Number(i2e = V a /u) 
because the gust function is r o /27rr 0 • V (a/r 0 ) in equation (3.86). But the acoustic 
pressure is more sensitive to the slope of the core 2/c(~ Re) shown in Figure 20, 
because again the behavior of the delta function, like the acoustic function, varies 
as the slope of the laminar core 2/c. 

In Figure 29, the nondimensional acoustic pressures for a point source, calcu- 
lated numerically based on the equation (3.151), are shown for two different values 
of turbulent vortex radius r 0 , nondimensjonalized by the half of the chord length 
6, with the fixed value of c and in the emission time domain s ( where s = 2 is the 
time required to pass through the chord). The different values give the different 
peak amplitude and width. Approximately, half of the value' of the vortex size 
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gives twice the value of the peak amplitude and half of the width, giving a sharper 
shape to the signal. This sharpened shape contributes to the higher frequencies of 
the acoustic pressure , as explained in Section 3.3, which are more annoying than 
the lower frequencies. So the quantity of interest here is the peak amplitude of 
the acoustic pressure, which is related to the shape of the signal. In Figure 30, 
the acoustic pressures are shown for the different values of c , inversely related to 
the Reynolds Number, with fixed vortex size r 0 . It is shown that the lower value 
of c gives the higher amplitude but the width is unchanged, which also makes the 
shape sharper. Considering that the pressure is nondimension alized by a quantity 
involving r o (~ 1/c), the pressure ( not nondimensionalized ) is more senstive to 
the value of c. 

For the peak acoustic pressure, the effects of the Reynolds Number and vortex 
size are shown in Figure 31. This shows that the peak acoustic pressure increases 
for the smaller vortex and the higher Reynolds Number. For this nondimension- 
lized peak acoustic pressure, it seems that the pressure becomes infite for the 
infinite Reynolds Number partially due to the property of the vortex and the be- 
havior of the point source in the incompressible flow. However, the pressure will 
remain finite for the noncompact source due to the interference effect between the 
sources. ? 

Noncompact Source and Observer Position Effects 

The results so far are for a point source and plotted in the emission time domain 
because it is not necessarily to introduce the time relation in equation (3.118) ( 
which makes the analysis simple) and the amplitudes are not changed between 


96 


the reception time domain, t , and the emission time domain , r , (as explained 
in Section 3.3 for a point source). The observer position is not considered in the 
results. However, for the noncompact source, the signal should be plotted in the 
reception time domain and the observer position should be considered as shown 
in equation (3.121) or (3.131). 

In Figure 32, the effect of the observer position on the acoustic pressure at the 
far field in the reception domain is shown for the spanwise noncompact source. 
The positions in the ij and i 2 directions are fixed as 5006 and the position in the 
x 3 direction ( spanwise direction) has the values 0, 1006 and 5006 for the spanwise 
source size of 206. In the figure, the effect of the position on amplitude and width 
of the signal for the fixed c and r 0 are shown. For the position at R 3 = 0, equation 
(3.121) is used because of the small interference between the sources and for the 
positions at J? 3 = 1006 and 5006, equation (3.131) is used. In both equations, the 
directivity function defined in (3.123) is not considered yet. So these effects in 
amplitude and width are mainly due to the interference of the noncompact source 
resulting from the integration in equation (3.121). 

In the case when the observer is in the plane of the symmetry of the source, 
that is f? 3 = 0, the peak amplitude is 17 in Figure 32 ; this value would be 20 
if there were no interference for the source size of 20 ( nondimesionalized by 6 ) 
because the peak amplitude of the point source for unit nondimensionalized source 
size is 1 in Figure 29 for c = 0.01 and r D = 4. In the case when the observer is far 
off the symmetry plane, for example i? 3 = 5006, the width of the signal becomes 
the order of the actual size of the source. Actually, the width in the reception time 
domain is almost the same as the time delay due to the source size, which is less 
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than the actual size of the source. This is explained in Section 3.3 ; the maximum 
time delay A t max (for a source Mach Number M = 0.5 and the same distance and 
source size) in Figure 22 is around 6, which is almost the same as the width of 
the signal, 6, in Figure 32. 

Directivity Pattern for the Peak Amplitude in the Acoustic 
Pressure 

To obtain the directivity pattern of the peak acoustic pressure, the observer po- 
sition is expressed by the angles 9,rp and the distance R* between the center of 
the source and the observer as defined in equation (3.146) and shown in Figure 
33. The distance R * is taken as 7146 corresponding to R x = 5006, R 2 = 5006 
and R 3 = 1006. Considering that the source is moving and the source strength 
, i.e., the unsteady lift, is a function of time, it is more reasonable to calculate 
the acoustic pressure for the moving observer with the same velocity as the source 
and having the same distance R * from the moving source in all directions. As 
shown in Figure 34, the interference effect ( the directivity function in (3.145) is 
not considered yet) due to the noncompact source is the same in the 9 direction 
in the plane of the symmetry for the moving observer because the time interval in 
reception domain is the same as the interval in emission domain in all directions, 
that is, = 1 in equation (2.33) for the moving source and moving observer with 
the same velocity. This time relation is one of the primary factors in the process of 
interference as explained in Section 3.3. On the other hand, the interference effect 
for the fixed observer depends on the 9 direction in Figure 34 because 
for M° = 0 in equation (2.33), that is, even though in the emission time domain 
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the time scales of the signals are the same, the time scales in the reception time do- 
main are different, depending on the location of the observer. However, as shown 
in the figure (not scaled), the effect of the moving observer is small for far field 
locations of the observer. 

The total directivity patterns for the peak amplitude including the directivity 
function defined in equation (3.145), in addition to the interference effect men- 
tioned so far are shown in Figure 34,35 and 36 for rp = 0,tt/ 6 and n/Z (xp is 
defined in Figure 33). As shown in Figure 34, at the plane of the symmetry of 
the source, the directivity pattern for the noncompact source is almost the same 
as that as the point source because of the small interference. However, for the 
plane at rp = ir/6 in Figure 33, there is a large difference in the directivity be- 
tween the noncompact source and the point source, as shown in Figure 35. This is 
mainly because the interference in the noncompact source in different 9 directions 
is different even though the time scales are the same since both the source and 
observer are in motion with the same velocity. At 9 = n/2 in the rp plane, the peak 
amplitude of the nonocompact source is the same as that for the point source, as 
in the case of rp = 0 because the interference is the same and is small. The pattern 
in rp = 7t/3 is almost the same as the pattern in ip — 7t/6 as shown in Figure 36. 

In summary, the directivity patterns of the peak acoustic amplitude are more 
concentrated in the plane of the symmetry compared to that of the point source. 
In that plane, the maximum peak acoustic pressure is observed and the value is 
almost the same as that of the point. source. 
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Noise Performance Trade-Off 

To reduce this maximum peak pressure , we need to reduce the vortex size r 0 and 
Reynolds number Re as shown in Figure 31. As explained in Section 3.4, one 
method would be to increase the distance, ( R — R m ) between the location of the 
maximum circulation and tip, because r 0 and Re both are directly related to it 
from equations (3.53) and (3.36). For given R, reducing Rm would result in a loss 
of lift near the tip which is proportial to (ClRm) 2 with the same lift coefficient. 
The other method would be to change the tip loading shape which is related to r 0 
for the fixed value of (R — Rm) and T 0 . 

Here, the influence of blade loading on the noise and performance has been 
studied for variations of the circulation shape of the form r/r„ = (sin0) 2m_1 
(m = 1 elliptic loading) as shown in Figure 37(b). For large values of m, the core 
radius, nondimensionalized by ( R — R m ) } increases and there is reduction of noise 
compared with the elliptic loading, but the tip loading relief also gives a reduction 
in the tip efficiency, defined in equation (3.183), as shown in Figure 37(c). By 
relating these two effects through the values of r 0 and c, Figure 38 is plotted. It 
shows that a maximum noise reduction of 5 dB in the peak amplitude ( which is 
related to the high frequency annoying noise) is obtained for m = 2 with 25% loss 
of tip efficiency and a noise reduction of 3 dB with 5% loss of efficiency. For the 
triangular loading there is more than 8 dB reduction in noise but also more than 
30% reduction in tip aerodynamic efficiency. 
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Chapter 5 
Conclusions 


The blade vortex interaction, where the turbulent viscous core is cut by a following 
blade, is analyzed under the assumptions that no distortion of the vortex path and 
no distortion of the vorticity of the core takes place during the interaction. From 
the analysis and results, several conclusions can be drawn as follows: 

(1) The detailed vortex structure, including the viscous and turbulent core 
generated by the rotor tip, can be approximately predicted with a simple turbulent 
core model and equivalent tip model. The size of the viscous core is of the order 
of the thickness of the blade and the size of the turbulent core is 3~5 times the 
blade thickness in the hovering case. 

(2) The interaction between a vortex core and a following blade is analyzed 
in the time domain , and a simple relation between the vortex structure and the 
unsteady pressure signal, using a gust analogy, is developed. 

(3) Acoustic pressure is obtained by introducing the acoustic response function 
for the step gust with a chordwise compact source. It shows that the spanwise 
non-compact source gives a sizable effect when the observer is off-center of the 
source because of the large interference effect. 
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(4) Maximum acoustic pressure varies with Reynolds Number( T 0 /v ) and 
vortex size, which is related to the tip loading distribution. The effect of a moving 
observer seems to be not significant in this blade vortex interaction. 

(5) Shape modification of the tip gives a maximum reduction of 6 — 8dB in- 
cluding the triangular loading ; a reduction of 3 dB in maximum acoustic pressure, 
compared with the elliptic loading, can be achieved for 5% loss in aerodynamic 
efficiency. 

For the extreme case of the rotor tip in which the blade passes through the 
vortex core, the chordwise non-compact source should be considered to obtain the 
correct high frequency effects. In the following remarks three limitations of the 
analysis are dicussed ; namely unsteady compressibility , noncompactness and the 
nonlinearity during the interaction. The vortex generator is also discussed. ( In 
the problem, three length scales are involved ; the wave length of the disturbance 
gust, the charateristic length of the body and the wave length of the acoustic 
signal observed.) 


r 
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Chapter 6 
Final Remarks 

Unsteady Compressiblity 

The exact linear unsteady compressibility effect is discussed by Graham^l for 
a purely sinusoidal gust. It shows that the effect, both in amplitude and phase, 
increases as the Mach number and the frequency increase. The linearized solution 
is still valid when the Mach Number becomes one if the frequency is sufficiently 
high^l , whereas the solution becomes singular for the low frequency as in the 
case of linear steady compressibility. Low frequency approximate solutions were 
obtained by Amiet^l, Kemp and Homicz^l. For the delt^ function gust, as 
shown by Amiet^l, the trend of the effect is similar to that for the purely sinu- 
soidal gust because the delta function has the same amplitude for low and high 
frequencies ( by definition). However, the effect of compressiblity is reduced for 
the step gust 1^1, compared to the above gusts, because the amplitude of the high 
frequency in the step gust is reduced. 

Thus, it is expected that the compressibility effect on the unsteady lift decreases 
more for the vortex gust, ( which is regular at the center of the vortex) because 
the amplitude of the high frequency reduces faster than that for the step gust. 
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The effect of compressibility for the point vortex passing beneath at the half 
chord height was shown by Adamczyk^l. and it would be useful to consider 
the above mentioned effects to determine the compressibility effect for the blade 
vortex interaction. 

Noncompactness 

The noncompact effect on the acoustic source appears when the characteristic 
length of the body is larger than the wave length of the acoustic wave radiated. For 
a purely sinusoidal gust, it seems that the effect of the noncompactness increases 
as the Mach Number and the frequency increase. However, assuming that the 
unsteadyness of the pressure near the leading edge dominates that near the trailing 
edge as used in high frequency analysis, the contribution to the total unsteady 
lift is confined at the small area near the leading edge. Then, the characteristic 
length to be used in the criteria for noncompactness should be the length of the 
small area in chordwise sense, which is much less than the geometric chord length 
usually used for the criteria. Therefore, it is expected that the critical frequency 
for noncompactness is reduced for the sinusoidal gust. 

For the vortex gust, the contribution of the high frequency is reduced as in the 
case of the unsteady lift and the contribution of the high frequency to the farfield 
acoustic pressure is reduced because the acoustic pressure varies inversely with 
the frequency. 

It is also expected that the effect is reduced more in the downstream direction 
because the observed frequency decreases due to the Doppler shift in frequency 
when the body is fixed in the wind ( as in the acoustic wind tunnel) or for a the 
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moving body with fixed observer. The effect is opposite for the upstream direction 
even though there is no difference between upstream and downstream for the body 
and observer moving with the same speed. The noncompactness in the upstream 
and downstream directions becomes more serious for the case of the interaction in 
cascade flow, and in the confined walls as in the case of fans and compressors of a 
turbo engine^, whereas it is less severe for the single vortex blade interaction. 

Although a clear criteria has not been given here for the effects of unsteady 
compressibility and noncompactness, the above comments for the blade vortex 
interaction are of interest. 

Nonlinearity 

The assumptions, used for the linear analysis, of the splitting of the velocity field 
into rotational and irrotational components, and the assumption of no distortion 
of the vorticity field or the vortex path, are valid for the case where maximum 
normal component of the gust velocity at the body surface is much less than the 
convecting velocity and where the thickness of the body is small compared to 
the wave length of the gust and large enough for no separation at the leading 
edge. In that case, the linear theory works well^l and it is known that splitting 
of the oncoming vortex by the body is the phenomena observed in experiments. 
For the nolinearity due to the thickness of an airfoil and the amplitude of a step 
gust velocity, Giesing^ showed that there is no significant effect on the Kiissner 
function for the 8% thick von Mises airfoil and for the gust velocity having 50% 
of a convecting velocity. 

However, for a free vortex interaction, the vortex path is distorted because of 
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the potential field of the vortex ; this is not to be expected for atmospheric gust a 
purely sinusoidal gust, a gust simulated by air injections and a weak blast wave. 
So, to predict the unsteady field quantitatively and to show the generation of the 
secondary vortex sometimes observed, the nonlinear and viscous effects should be 
considered. The quadrupole effect due to the nonlinearity in the flow field, which 
would be the same order as the dipole, and the transonic effect, vortex induced or 
not should be also considered in the case of higher compressible flow. Even in the 
incompressible case, the development of the prediction scheme for high Reynolds 
Number flow with reasonable vortex core is desirable. 

The reasons for the large variation at the leading edge of the blade, observed in 
the data from the large scale experiments both at low and high Mach Number^ 93 ', 
are still not revealed. It might be due to the high frequency effect discussed in 
the high speed linear theory or the forward propagation of the shocki 3 ®! or the 
distortion of the vortex and vortex path. 

Vortex Generator 

In addition to those for the interaction, it is of interest to comment on the vortex 
generator since the vortex strength and structure are important factors for the 
interaction. In the measurement of a vortex flow, it is very useful to measure 
both quantities for the vortex and the vortex generator whether it is fixed or 
rotating (,t.e. maximum circulation and its position, circulation distribution and 
aspect ratio of the generator) because they are important parameters to predict 
the vortex profile and to determine the distance of fully roll up vortex before the 
interaction ( the maximum circulation distance is changed for different tip shapes 
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and tip speeds in the case of the rotating blade). It should be also pointed out 
that for a two dimensional vortex, generated by an oscillating airfoil or an airfoil 
in impulsive motion, that the behavior is different from the tip vortex because of 
no axial velocity in the 2-D vortex. However, it might simulate the tip vortex in 
the decay stage where the axial velocity is negligible. 

Care should be taken in comparing the vortex theory, which usually assumed 
full roll up of the core, with the experiment data because sometimes it is observed 
that the maximum circulation of the concentrated vortex is half of the circulation 
at the root of the generator^. 
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Appendix A 


A.l 


The Heaviside function has the following properties 


mn 

df 


= 6{f) 


(A.l) 


± i{ H( f )F } = H( f) ^ + Fi mi (A. 2) 

where F is an arbitrary function and 8(f) is the delta function having zero value 
except at / = 0, that is, at the surface. 

Using the equation of the surface motion and the above properties, the follow- 
ings can be obtained 

dH(f) = dH(f) df 
dt dfdt 

= S(/) ■(-«§£) (A.3) 
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and 


dHjJ) 

dxi 


dH df 
df dx{ 

«/>■•£ 


(A.4) 

(A.5) 


A. 2 


The volume integration having the delta function in the integrand can be changed 
into the surface integration 

= (A6) 
where S is the surface of the body in the motion defined as 

/(*?., t) = 0 (A.7) 

The time integration having the delta function, as shown below, in the integrand 
can be carried out as 


j'6{g{T))F{T)dT = 


nil 

dg{ri)/dT | 


where r e ‘ is the i th root of 


ff(K) = o 


(A- 8) 


(A- 9) 
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Appendix B 


The equation is 

— + G , '+4p 2 G=C 1 (B. 1) 

z 

with the boundary condition 

G - 0, z = 0 (B. 2) 

An appriximate solution in the following form is tried 

G = c ■ g 0 {r)) + c 2 ■ gx(ri) (B.3) 

where rj is a stretched coordinate given by 

V = * 2 /c 2 (BA) 


Substitution into the equation gives equations for go(f}) and gi(r;), i.e. 

4i]g' 0 + 2 < 7 q — Ci = 0 (B-5) 
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From equation (B.5) 


9[ + 4ri9o9i + 1P 2 9 o - 0 


(B. 6) 


</TT4C^j-l 

9o — A _ 

4r) 


(B. 7) 


After integration, the solution for g 0 satisfying the boundary condition is obtained 


as 


VI + 4 Ctf 1 1 , (\ + VI + 4^17 

So = 2 2 _ 2 ln i 


(B.S) 


From equation(B.6) 


, -'Lp 2 9o 

9i = 


1 + 4r)9o 


(B.9) 


Using (B.8) and (B.7), the above equation becomes 


9l = -P 2 1 - 


1 In (khZHl 35) 

s/T+^C^rj y/T+~4crf 


(B. 10) 


, which can be integrated to yield 


/? 2 ( 1 + Vl + /" 1 + \/l + 4C : 

= J '“I 


S)j 


(5-11) 


ill 


Appendix C 


L(s) = 2icbpU b 


2nr o {0) 


v&m+J' 7 vi 

T n J 0 T q 



ip(s — (j)dcr 


where 

x{j(s) = 1 - 0.5e" 013 ' - 0.5e~* 
rp(s — cr) = 1 - O.5e _013, e 013<r - 0.5e“V 


I(s) = L(s)- 


1 27rr o (0) 

2icbpU b T„ 


0 < s < Si 


I(s) 


= -~rp{s)+ f -- — br ° rp(s - cr)dcr 

x 0 Jo ( x 0 - crb) 1 2 


**.y 


_ro_ _ ro r 1 

x 0 b Jo {cr - 

+ + 0.5 y f 

x 0 b Jo (a — -f) 2 

= Pi{s) + P 2 {s)e~ 013 ' + P z {s)e- 


-- + o.5 r f f i 

x Q b Jo (cr — 


,0.13<t 


£s02 
6 ) 


da\ 


0.13* 
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where 
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0 {* - t ) 2 


e ko da = 


-1 


„ka 






[• e ka 

Jo a — 


-da 
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a ko 


k*L f k W •) e ' x a 

e » / ——dt , £ = — A:(cr — — ) 

Jk?f t 'b 

[r ^j-dt~r~dt\ 

[ * J k^ t j 

{*(**-.))-*(**)} 


o a 


* a. 
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k £ 2 . 

= — e * 


Jfc £ 2 . 

= — e * 


Sl ^ 3 ^ S2 


/(*) = -%*(•)+ J t 
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(x 0 - ab ) 2 

Pi( s i) + P 2 (sj)+ 0 . 5 — f 

r» rn J t\ ( a 


/■» b 

rp(s-a)da+ / ip(s - a)da 

J *i r 0 
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r o J‘> (<r-f ) 2 

p ,M + 0 


s 0A5tT da 


,-0.13* 


= P 1 ( S )+P 2 ( S )e- 013 ' + P 3 ( 5 )e-' 


Si < s < s 2 


I(s) = — — V'(-s) -J- / 

Xq Jo 


br n 


(x 0 - ab ) 2 


ip(s-a)da+[ - — ip(s-a)da 
J>i r„ 


Pl(Sl) ~ 7 0 \ + [ p ^ + 0 5 7 0 L {^W e °' 


_ b 
»1 >0 

QAZa da 


,-0.13* 


= PM + P 2 (s)e- 01U + P 3 (s)e- 
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S 2 < S < S 3 


= -^ (s > + r -(^bf p( ~ s - ° )ia + r -£*<* - ^ 

+ / -^{s-<j)d<j 

J 82 C T 0 

-- |pi(s 2 ) + ~] + |p 2 (s 2 ) - 0.5~ J' e 01z °da e" 013 ' 

2 b t* 

+ Ps(s 2 ) — 0.5 / e^da e * 

C fg 


= Pi(5) + P 2 (s)c-° 13 * + P 3 (s)e-* 
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where 


where 


+ [*(-> + .o • 5 T/..'(^ e ^ 

Pi(s)+P 2 (s)e- 013 * + P 3 (s)e-' 


i: 


a ka 


■da = e k *? 


¥ £^ ( , = - k(a - % 
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Appendix D 
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5j ^ 3 ^ 32 


dPijs) _ _b_ dP 2 {s ) = dP 3 (s ) 
da r„ ’ ds ds 

s 2 < s < s 3 

gPt(s) _ 2 b dP 2 (s) = dP a (s) 

ds cr 0 ’ ds ds 
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A: = 1 ; » = j = 2, k = 0.13 ; i = j = 3 
Table 1: 
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I developing stage 

II developed or persistence stage 

III decay stage 

IV unstable stage 


Figure 2: Stages in develoment of tip vortex 





No microphone 
response 



Figure 4: Spectrum for helicopter noise (UH-1A, Tiedown Thrust 6001b, Tip 
Velocity 720ft/sec, Microphone distance from source 200ft) 
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Figure 5: Impulsive noise signal 
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Figure 6: Spectrum of rotational noise^l (a) with blade vortex interaction (b) 
with tip vortex eliminated (UH-1H, 1.5 g left turn, (i = 0.24, 120mph ) 
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Figure 7: Trajectories of rotating tip vortex^' (a) side view (b) plan view 
(j, = 0.164 
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Figure 13: Vortex velocity distribution 
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Figure 14: Coordinate system of the blade during interaction 
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Figure 15: Kiissner function 
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Figure 16: Gradient of the Kiissner function 
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Figure 21: Interference process of the spanwise noncompact source 
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Figure 22: Variation of delayed signal time with observer position 



Figure 23: Simplified model signal in time domain and its spectrum 
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Figure 25: Experimental data for rotor tip vortexl 42 ! 
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Figure 26: Experimental data for rotor tip vortext 4 ^ 


147 




Figure 27: Vortex core velocity profile comparition (a) circulation profile (b) tan 
gential velocity profile 
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Figure 29: Variations of signal with vortex size 
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Figure 30: Variations of signal with Reynolds Number ( T 0 /v ) 
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Figure 32: Variations of signal with observer position 
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Directivity Pattern of Peak Acoustic Pressure 
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Figure 34: Peak acoustic pressure in the plane of symmmetry 
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Figure 35: Peak acoustic pressure in the plane rp = W6 
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Figure 36: Peak acoustic pressure in the plane rp = tt/3 
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Figure 37: (a) Tip circulation profile (b) Variation of r 0 /R — R m and A p(dB) with 
tip- loading parameter m (c) Variation of tip efficiency with tip-loading parameter 
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Figure 38: Trade off between noise reduction and aerodynamic efficiency for vari 
ous m. 
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